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This  report  was  prepared  by  the  General  Applied 
Science  Laboratories  (GASL)  under  the  terms  of 
Contract  DAAJ02-68-C-0092 .  It  consists  of  a 
discussion  and  presentation  of  findings  along  with 
the  approach  followed  to  design,  fabricate,  and  test 
a  high-speed  burner  for  turbomachine  applications. 

The  unique  feature  of  a  high-speed  burner  is  that  the 
combustion  process  can  take  place  at  higher  Mach 
numbers  than  it  can  with  current  gas  turbine  combustion 
systems.  Being  able  to  conduct  the  burning  process 
at  higher  Mach  numbers  will  reduce  the  diffusion 
requirements  of  a  compressor  diffuser,  therefore 
resulting  in  a  more  efficient  compressor  diffuser 
and  combustor  system. 

The  object  of  this  contractual  effort  was  to  determine 
the  feasibility  of  a  high-speed  combustion  chamber 
having  a  burner  entrance  Mach  number  of  approximately 

0.5. 

In  general,  this  experimental  effort  resulted  in  the 
demonstration  of  ignition  and  partial  combustion  of 
liquid  JP-4  fuel  with  burner  entrance  Mach  numbers 
as  high  as  1.0. 

The  conclusion  and  recommendations  contained  herein 
are  concurred  in  by  this  Command.  This  concurrence 
does  not  imply  that  a  high-speed  burner  is  practical 
for  gas  turbine  applications.  However,  it  is  believed 
that  the  concept  shows  some  promise,  and  further 
investigations  will  be  required  to  prove  or  disprove 
the  practicability  of  such  a  system. 
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ABSTRACT 


The  theoretical  analyses,  the  design, and  the  experimental 
verification  of  a  high-speed  combustion  chamber  are  described. 
For  turbomachines ,  this  type  of  burner  is  used  when  compressor 
outflow  speed  is  so  high  that  diffusion  to  low  speed  presents 
severe  pressure  loss  penalties.  The  present  study  showed 
that  for  a  low-mass-flow,  high-pressure-ratio  turbomachine, 
combined  diffuser  and  combustor  losses  are  minimum  for  a 
burner  entrance  Mach  number  of  about  0.5.  The  use  of  the  GASL 
finite  rate  chemistry  and  turbulent  mixing  programs  is 
discussed  along  with  the  combustor  modeling  and  flame  spread 
predictions.  Finally,  a  series  of  experiments  is  described, 
and  burner  pressure  loss  and  temperature  profiles  are  shown 
over  a  wide  range  of  burner  airflow  conditions,  i.e.,  pressures 
from  1  to  11  atmospheres  and  inlet  air  temperatures  from 
ambient  to  1200°R. 
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I.  INTRODUCTION 


There  are  several  current  trends  in  modern  gas  turbine  design, 
particularly  in  the  small  power  range,  that  are  creating 
difficulties  for  the  combustion  chamber  designer.  For  in¬ 
stance,  the  trend  toward  higher  cycle  temperatures  not  only 
increases  the  required  overall  equivalence  ratio  for  the 
burner  but  also  increases  the  specific  power  of  the  machine. 

The  resulting  smaller  dimensions  of  the  rotating  parts  tend 
to  make  the  combustion  chamber  dominate  the  engine  packaging 
envelope  when  the  conventional  recirculating  concept  is  used. 

In  addition,  the  increased  cycle  temperature  requires  an  in¬ 
creased  cycle  pressure  to  lower  the  specific  fuel  consumption. 
The  foregoing  considerations  combined  with  the  power  re¬ 
quirements  for  Army  applications  constrain  this  engine  con¬ 
figuration  to  include  a  centrifugal- type  compressor.  The 
outlet  flow  from  such  a  compressor  wheel  has  very  high  velocity; 
at  the  pressure  ratios  under  consideration  (10  to  12),  the 
flow  is  supersonic.  Thus  the  diffusing  stator  in  this  appli¬ 
cation  is  an  extremely  critical  element,  in  particular 
since  diffuser  performance  drops  off  drastically  at  higher 
Mach  numbers. 

One  way  to  improve  performance  is  to  consider  the  diffuser 
and  combustor  as  a  system  and  to  trade  off  reduced  diffusion 
in  the  stator  for  increased  burner  inlet  airspeed.  However, 
the  efficiency  of  the  conventional  recirculating  primary-type 
burner  is  unacceptably  low  at  airspeeds  much  above  200  ft/sec, 
so  that  a  fundamental  change  in  burner  design  philosophy  is 
required  for  such  a  system. 

The  purpose  of  the  present  investigation  is  to  utilize  the 
philosophy  and  some  of  the  computational  technology  developed 
for  SCRAMJET  (Super sonic-Conbustion-Ram-Jet)  engines  to  design 
this  new  type  of  turbo-engine  burner.  Accordingly , a  GASL 
computer  program,  which  analyzes  flows  with  finite-rate 
combustion  of  hydrocarbon-type  fuels  plus  turbulent  mixing, 
was  utilized  to  determine  the  design  layout.  No  particular 
engine  was  specified,  since  this  work  was  of  fundamental 
nature;  indeed, one  of  the  tasks  of  the  program  was  to  determine 
the  inlet  conditions  for  best  performance.  However,  some 
guiding  parameters  were  set  initially  by  the  Army  on  mass 
flow  range,  pressure  ratio,  turbine  inlet  temperature, and 
overall  size,  so  that  the  results  would  be  applicable  to  their 
future  engine  needs.  These  engine  parameters  are  listed  in 
Table  I. 
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TABLE  I.  ENGINE 

PARAMETERS 

Sea  Level 

Sea  Level 

Military  Rated  Power 

Idle* 

Engine  Airflow  Rate,  lb/sec 

5.00 

2.96 

Compressor  Pressure  Ratio 

12.5:1 

5.08:1 

Compressor  Adiabatic  Efficiency;  % 

80 

78 

Turbine  Inlet  Temperature  ,°F 

2500 

1587 

Maximum  Fuel  Inlet  Temperature, °F 

300 

300 

Absolute  Altitude  Limit, 25,000  ft 

Ram  Pressure  Ratio  =  1.0 

*  20%  Power 

Although  it  was  realized  that  the  burner  had  to  perform 
efficiently  over  the  entire  performance  spectrum  indicated  in 
Table  I,  it  was  decided  that  the  sea  level  100  percent  power 
condition  would  be  most  important.  This  case  was  called  the 
design  point,  and  all  the  major  design  calculations  were  per¬ 
formed  at  this  condition. 


The  report  is  divided  into  sections  as  follows:  Section  II 
gives  the  analysis  which  determined  the  optimum  burner  inlet  ^ 
Mach  number  at  the  design  point  and  thus  fixed  the  inlet  and 
exit  flow  areas.  The  performance  variation  is  then  computed 
over  a  range  of  power  settings  both  at  sea  level  and  at 
altitude.  Section  III  describes  the  burner  design,  which 
includes  fuel  injection,  ignition,  and  the  finite  rate  chemical 
reaction  system.  Section  IV  points  out  salient  features  of 
the  experimental  facility,  while  Section  V  presents  and 
discusses  the  experimental  results.  Conclusions  and  rec¬ 
ommendations  are  given  in  Section  VI. 


II.  BURNER  INLET  CONDITIONS 

In  order  to  determine  the  optimum  burner-inlet  Mach  number, 
the  type  and  performance  of  the  burner  and  also  of  the  diffuser 
must  be  known.  In  this  section,  the  performance  of  three 
candidate  burner  configurations  is  reviewed  on  the  basis  of 
one-dimensional,  equilibrium  flow.  A  dump-diffusion  process 
is  chosen  for  the  diffusion  model,  and  the  overall  diffuser- 
combustor  performance  is  computed  over  a  range  of  inlet  Mach 
numbers.  The  optimum  burner  inlet  Mach  number  is  then  defined 
when  the  overall  pressure  loss  is  minimum. 

Before  showing  these  analyses,  definition  of  the  off-design 
turbomachine  flow  is  given  so  that  physical  flow  quantities 
can  be  associated  with  each  power  setting. 

OFF- DESIGN 

The  following  assumptions  are  used  in  defining  the  off-design 
flow  conditions: 

1.  Compressor  operating  line  is  synthesized  from  the 
"back -bone"  type  correlation  detailed  in  Reference  1. 

2.  Combined  diffuser-burner  stagnation  pressure 
loss  equals  0.81. 

3.  The  flow  is  choked  in  the  first  turbine  nozzle 
stage. 

4.  Turbine  adiabatic  efficiency  (both  gas  generator 
and  power)  equals  0.85. 

5.  Complete  expansion  occurs  in  the  power  turbine. 

Table  II  lists  the  mass  flow,  burner  and  turbine  inlet  stag¬ 
nation  temperatures,  and  diffuser  inlet  (rotor  discharge) 
stagnation  pressure  over  a  range  of  output  power  settings. 
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TABLE 

II.  OFF- 

-DESIGN  CONDITIONS 

Percent 

Altitude 

Power 

Sea  Level 

25.000  ft 

Mass  Flow,  m 
(lbm/sec) 

20 

60 

2.96 

4.37 

1.20 

1.65 

100 

5.00 

2.0 

Burner  Inlet 

20 

970 

795 

Temperature  T 

60 

1120 

930 

(  °R  ) 

100 

1190 

990 

Turbine  Inlet 

20 

2046 

1690 

Temperature  T 

60 

2640 

2150 

(  °R  )  4 

100 

2960 

2440 

Diffuser  Inlet 

20 

94 

35 

Total  Pressure,  P2 

60 

150 

52 

(psia) 

100 

190 

70 

These  cases  are  studied  in  the  following  section  on  the 
estimation  of  combustor  performance. 

ESTIMATION  OF  COMBUSTOR  PERFORMANCE 

The  performance  of  three  burner  models  is  examined  in  this 
section  to  determine  the  performance  of  each,  i.e.,  pressure 
recovery  and  outlet  to  inlet  area  ratio  in  terms  of  the  inlet 
air  Mach  number.  The  models  chosen  were  constant  pressure, 
constant  Mach  number/  and  constant  area  burning.  Equilibrium 
chemistry  was  assumed  along  with  quasi  one-dimensional  flow, 
and  no  account  was  made  of  the  fuel  addition  processes. 

Constant  Pressure  Combustor 

Conditions  at  the  end  of  burning  are  determined  for  a  con¬ 
stant  pressure  burning  process.  The  calculations  are  made 
assuming  a  constant  molecular  weight  and  values  of  y  -  1.4 
at  the  start  of  burning  and  y  =  1.3  at  the  end  of  burning. 
Experience  at  GASL  has  shown  that  these  assumptions  give 
very  good  results  for  an  equilibrium  combustion  process  in 
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the  range  of  fuel-air  ratio  of  interest  in  the 
cation. 

present  appli 

Constant 

pressure  burning  is  defined  by 

and 

P3  -  P4 

(1) 

where 

< 

u> 

II 

(2) 

p  =  static  pressure 

V  =  air  velocity 
subscript  3  =  before  burning 
subscript  4  =  after  burning 

From  the 

second  condition, 

(V/a* ) 4  =  (V/a*)3  a*3/a*4 

(3) 

where  a*  is  the  critical  speed  and  (V/a*)  is  a 
Mach  number  given  by 

function  of 

V/a*  =  M  /  y  -+-  *  ~ 

V2  (1  +  M2) 

(4) 

where 

M  =Mach  number 
y  =ratio  of  specific  heat 

Now 

a*  =  V>'gRT* 

(5) 

and 

2 

T*  =  l  +  y  T 

(6) 

so  that 

a*  *  Vrfry  gRT 

(7) 

where 

T*=  critical  temperature 

T  =  total  temperature 
g  =  acceleration  of  gravity 

R  =  gas  constant  =  1544/mw 

Hence,  Equation  (3)  may  be  written 


i 


kL 
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(8) 


(V/a* ) .  =  (V/a*) 


y3(Uy4)T3 
3  '  y4(1+y3,T4 


For  given  values  of  and  T3,  this  equation  gives  the  varia¬ 
tion  of  M3  with  M4  as  shown  in  Figure  1.  Since  the  ratio 
T3/T4  is  more  or  less  independent  of  altitude,  all  of  the 
present  results (except  the  amount  of  heat  added)  are  only 
functions  of  power  setting. 

Equations  (2)  and  (3)  applied  to  the  continuity  equation  give 


where 


A  =  area 

t  = static  temperature 


This  may  be  written 

A4  ^  (t/T), 

A3  =  T3  (t/T) 3 

where  t/T  is  the  Mach  number  function  given  by 


t/T  =  1/(1 


y-1 

+  m 


It  is  clear  that  with  M4  determined  from  Equation  (8)  , 
Equation  (10)  gives  the  burner  area  ratios  shown  in  Figure  2. 

The  burner  recovery  is  defined  by 


= 

where  P  is  total  pressure. 

Since  P4  =  P3,  this  may  be  written 

(p/P)3 

VB  =  (p/P)4 

where  p/P  is  the  Mach  number  function  given  by 
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p/p 


1 


Izl 

2 


_L_ 

y~  i 


(13) 


With  M4  determined  from  Equation  (8),  Equation  (12)  gives 
the  burner  recoveries  shown  in  Figure  3. 


The  determination  of  the  curves  of  fuel-air  ratio  in  Figure  4 
was  made  using  the  data  in  Reference  3.  These  data  are  for 
the  general  hydrocarbon  CnH2n#  but  the  results  apply  to  JP-4 
as  well. 


Constant  Mach  Number  Burner 


The  constant  Mach  number  burning  process  is  treated  by  con¬ 
sidering  the  general  relationship  between  total  pressure  and 
total  temperature,  which  for  a  one-dimensional  process  is 
expressed  by 


dp  yM2  dT 

P  2  T 


(14) 


In  the  present  application,  M  is  constant  but  y  varies.  How¬ 
ever,  since  the  change  in  y  is  small,  it  is  reasonable  to 
treat  y  as  constant  at  <.  i  average  y: 


y 


y3  +  y4 
2 


Then  Equation  (14)  gives 


(15) 


(16) 


The  static  pressure  ratio  across  the  burner  is  determined 
from 


P4  <P/p)4 

p3  '  (p/P)3 


(17) 


(Note  that  although  M4  is  equal  to  M3,  (p/P)4  is  not  equal  to 

(p/P) 3  because  of  the  difference  between  y ^  and  y3 . ) 


The  temperature  ratio  and  velocity  ratio  are  given  by  the 
general  relationships 
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Altitude,  Thousands  of  Feet 


Figure  4. 


Fuel-Air  Ratio 
Power  Levels. 


vs  Altitude  at  Various 
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3 
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( 


_2 _ 

V1 


2 


)  M3  ) 

)  mJ) 


The  area  ratio  may  then  be  found  from 


V  t 

3  4 

>  <^r>  <  r> 

4  3 


(18) 


(19) 


(20) 


For  the  constant-Mach-number  burning  process,  the  results 
presented  are  the  burner  recovery  (Figure  5) ,  the  area  ratio 
(Figure  6),  and  the  pressure  ratio  (Figure  7)  . 


Constant-Area  Combustor 


The  constant- area  burning  process  is  defined  by  the  conserva¬ 
tion  of  total  momentum  as 


where 


p  A  +  —  V.  =  p.A  +  —  V, 
^4  g  4  *3  g  - 


A  =  area 
m  =mass  flow 

g  = acceleration  of  gravity 


(21) 


From  the  continuity  conditions, 


pA  = 


mRt 

V 


(22) 


and  the  molecular  weight  is  assumed  to  be  the  same  before 
and  after  burning. 

Using  Equation  (22)  in  Equation  (21)^ 


mRt. 


V. 


S  v  = 

g  4 


mRt_ 

- i  +  5V 

v,  g  3 


(23) 
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Burner  Recove 


Figure  5.  Burner  Recovery  -  Conrtant-Mach-Nunvber 
Burning. 
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Burner  Area  Ratio,  A4/A3 


Figure  7.  Burner  Pressure  Ratio  -  Constant-Mach-Number 
Burning. 
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or 


Using  Equation  (11),  the  velocity  is  written 


V  = 


(25) 


From  the  definition  of  Mach  number, 


_  _1_ 

2  2 
V  yM 


(26) 


and  Equations  (24)  and  (25).  the  equation  defining  constant- 
area  burning  is  written 


where 


F  = 


2  2 
M  )  M 

~2 


+  yM 


X 


(27) 


(28) 


It  is  apparent  that  with  M-,  T  and  T.  known.  Equation  (27) 
gives  m4<  j  j  4 


With  M4  determined,  the  pressure  after  burning  is  found  from 
the  condition  of  conservation  of  total  momentum. which  is  now 
written 


or 


P4V1+V4  )  =  P3V1+V3) 

P4  1  +  y3«3 


(29) 

(30) 
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The  burner  recovery  (ratio  of  total  pressure  after  and  before 
burning)  is  determined  from 


V 


B 


iiU 

<  f  >4 


(31) 


where  (p/P)  is  the  Mach  number  function.  Equation  (13). 


The  results  of  calculations  for  the  constant-area  burning 
process  are  given  as  the  burner  recovery  (Figure  8),  the 
Mach  number  after  burning  (Figure  9),  and  the  pressure  ratio 
(Figure  10) . 


Comparison  of  the  constant. Mach-number  and  constant-area 
burner  with  those  for  the  constant- pressure  burning  indicates 
the  superiority  of  the  constant-pressure  burning  process; 
namely,  the  burner  recovery  is  higher.  A  quantitative  compari¬ 
son  of  the  burner  recovery/ at  100  percent  power,  is  shown  in 
Figure  11.  The  results  also  show  that  the  higher  performance 
requires  a  larger  burner  exit  area  and  results  in  lower  burner 
exit  Mach  number. 

DIFFUSER  -  TURNING  LOSS  MODEL  AND  EVALUATION  OF  OPTIMUM 
BURNER  ENTRANCE  CONDITIONS 


The  optimum  burner  entrance  conditions  are  defined  when  the 
pressure  recovery  from  compressor  outlet  to  turbine  inlet  is 
maximized.  In  order  to  calculate  this,  the  efficiencies  of 
each  element  included  in  the  system  must  be  determined.  For 
the  combustor,  the  pressure  recovery  has  already  been  defined 
as  a  function  of  inlet  (to  the  burner)  Mach  number  in  the 
constant-pressure  burner.  The  associated  loss  picture  for  the 
diffuser  model  and  turning  duct,  since  we  are  considering  a 
radial  compressor,  is  now  required. 

This  section  describes  such  a  model  and  its  resulting  pres¬ 
sure  recovery.  From  the  combination  of  burning  and  diffuser¬ 
turning  pressure  loss  as  a  function  of  burner  inlet  Mach 
number,  the  optimum  operating  condition  is  determined. 

Diffuser  Model 


There  are  two  extreme  examples  of  a  diffuser  model:  first, 
it  can  be  perfectly  efficient  so  that  the  optimum  burner 
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Burner  Exit  Mach  Number 


*  ”  .  -J 

Burner  Entrance  Mach  Number,  m3 


Figure  9.  Burner  Exit  Mach  Number  -  Constant-Area  Burning. 
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Burner  Pressure  Ra 


Figure  10. 


Burner  Static  Pressure  Ratio 
Burning. 


Constant-Area 


inlet  Mach  number  is  zero;  or  second,  the  diffuser  can  be 
completely  inefficient  (so  that  all  the  kinetic  energy  in 
the  flow  from  the  compressor  is  lost)  and  the  optimum  burner 
Mach  number  equals  the  compressor  outlet  value.  Of  course, 
the  actual  diffuser  is  neither  of  these,  but  the  prediction  of 
its  performance  is  still  one  of  the  most  difficult  fluid 
dynamic  problems  yet  to  be  resolved. 


Nevertheless,  in  order  to  analyze  a  simple  but  representative 
diffuser  system,  the  following  items  are  considered:  (1)  the 
Army's  requirement  of  low  weight  and  volume;  (2)  the  fact  that 
separated  flow  exists  in  most  diffusers;  and  (3)  any  simple 
model  should  be  length  or  angle  independent.  These  items 
lead  one  to  choose  a  dump-diffuser  as  shown  schematically  in 
the  sketch  below: 


Station  2/3 


Station  2a 


Rectangular  Diffuser  and  Elbow 

Rectangular  piping  is  indicated  because  it  is  representative 
of  a  centrifugal  compressor  application. 

To  define  conditions  at  station  (2a)(i.e.,  just  upstream  of 
the  dump  diffuser),  we  must  describe  the  process  whereby  the 
supersonic  outflow  from  the  compressor  is  diffused  below 
Mach  1.  The  most  efficient  actual  diffusion  is  to  assume 
that  the  flow  experiences  a  normal  shock.  This  would  corres¬ 
pond  to  a  long  pipe  frictional  type  diffuser.  From  the  normal 
shock  relations,  Cf,  Reference  4,  for  M2  =  1.27,  we  find 
M2a  =  0.8  and  P2q/p2  =  °*987.  From  dump  loss  data  shown  in 
Reference  5,  conditions  at  station  (20)  are  defined.  The 
turning  losses  are  determined  by  using  the  data  shown  in 
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Reference  6  and  assuming  no  area  change  in  the  turning  duct. 
Thus, the  conditions  at  the  burner  inlet  are  defined. 

Analysis 


To  find  the  dump-diffusion  pressure  loss  from  Reference  5 , 

_  _  r  _  — «  2 


where 


'  2a  R 

(32) 

^2  a  “U  J 

1  2 

2  ^  P  ,  the  dynamic  head 

(33) 

=  A.  /A_„ ,  the  area  ratio 

2/3  2a 

(34) 

and  z  is  a  compressibility  factor. 


Equation  (32)  can  be  rewritten,  substituting  Equations  (33) 
and  (13 )  as 


=  i-  \  ^ 
P2tt  2  2“ 


l+( 


y- 


i  r a -i  “| 

z  -J— 

-  ar  - 


(35) 


and  for  M2  =  0.8 


^  =  x  .  o. 


304  z  ( 


V1  .2 


(36) 
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The  value  of  z  for  a  rectangular  diffuser  at  this  Mach  number 
(0.8)  is  2.6. 


To  find  the  Mach  number  after  the  dump  (at  station  (2/3)  ), 
we  note,  since  the  stagnation  temperature  of  the  air  is 
constant,  that 


*V  )( 

°  2/3  R 


p/P 

^  iJt/T  ' 

^(t/T)V 

2a 

p/P 

. 

2/3 


(37) 
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where  t/T  and  p/p  are  defined  by  the  Mach  number  functions. 
Equations  (11)  and  (13)  respectively. 

Comparison  of  the  pressure  recovery  P2$/p2a  f°r  t*ie  ^umP 
diffuser  with  some  diffuser  experiments,  Reference  7,  is 
shown  in  Figure  12.  The  experimental  data  are  for  two 
length-to-width  ratios  (L/W)  diffusing  from  Mach  1  to  M2/j* 
These  data  show  the  rather  typical  increase  and  then  decrease 
which  is  due  to  flow  separation  and  skin  friction.  As  L/W 
decreases  toward  0  (a  dump  diffuser) ,  it  is  seen  that  the 
experimental  curves  trend  toward  the  dump  value. 

To  determine  the  pressure  loss  in  the  turn,  from  Reference  6, 


or 


P2/3~  P3 
'  q2/* 


=  c 


(38) 


1-t 


y  2 

M-o 
2  2p 
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y-i 


(39) 


For  a  rectangular  duct  whose  inlet  perimeter  has  a  width 
twice  its  depth  and  a  radius  of  curvature  equal  to  the 
width,  the  value  of  C  =  0.2. 


The  pressure  loss  for  the  entire  diffuser-combustor  system 
is 


V 


<  r  >  (V 

v2  v2a  2/3 


(40) 


where  T)  is  the  combustion  chamber  pressure  loss  shown  in 
Figure  3  as  a  function  of  Mach  number  and  power  setting.  By 
denoting  M3  as  the  design  Mach  number  defined  as  the  sea  level 
100-percent  military  power  condition; the  value  of  r?B  is 
determined.  A  plot  of  Equation  (40)  as  a  function  of  burner 
inlet  Mach  number  is  shown  in  Figure  13,  where  it  is  seen  that 
X)  first  increases  and  then  decreases  with  increasing  M3.  There¬ 
fore,  the  Mach  number  where  77  is  maximum  is  the  optimum 
burner  inlet  condition  and  is  approximately  0.5. 


The  dotted  portion  of  the  curve  shows  the  recovery  when  the 
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Diffuser  Recovery  -  For  an  Expansion 


experimental  diffuser  (Reference  7,  L/W  =  10)  values  are 
used  instead  of  the  dump  diffuser  values.  The  recovery  of 
the  system  in  this  range  of  experimental  data,  although  greater 
than  with  the  dump  diffuser,  is  still  less  than  the  maximum 
recovery  with  the  dump.  This  means  that  if  experimental  data 
were  available  for  this  high  value  of  M 2 it  would  show  that 
the  optimum  Mach  number  is  between  0.53  and  0.33  (the  better 
the  diffuser,  the  lower  the  burner  inlet  Mach  number)  .  There¬ 
fore,  it  still  appears  reasonable  to  choose  the  design  Mach 
number  equal  to  0.5. 

METHOD  OF  COMPUTING  BURNER  OFF-DESIGN  PERFORMANCE 

The  combustor  loss  curves  are  essentially  design  point  values. 
That  is,  a  given  combustor  inlet  Mach  number  not  only  defines 
the  burner  outlet  to  inlet  area  ratio  (A4/A3)  but  also 
specifies  li.c  combustor  inlet  conditions,  i.e.,  area,  velocity, 
static  temperature,  etc. 

To  determine  the  performance  of  the  design  Mach  number 
combustor  over  the  power  spectrum,  the  following  model  for  the 
off-design  combustion  process  is  postulated: 

1.  Burning  occurs  at  constant  pressures  at  all 
power  settings.  The  associated  burner  recovery 
and  requisite  flow  area  ratio  are  thus  defined. 

2.  Burning  is  initiated  at  the  combustor  inlet 
station. 

3.  isentropic  compression  occurs  from  the  end  of 
burning  to  the  actual  exit  station. 

The  reasons  for  choosing  such  a  burning  model  are  as  follows: 

1.  It  is  fairly  representative  of  the  physical 
situation  in  the  combustor. 

2.  Pressure  losses  can  be  defined  easily. 

3.  Determined  pressure  loss  is  conservative. 

To  clarify  this  last  point,  if  combustion  would  take  place 
within  the  entire  chamber,  then  the  average  Mach  number  at 


I 


which  burning  occurs  would  be  less  than  the  one  in  the  present 
model.  Since  losses  are  directly  proportional  to  Mach  number, 
the  burning  losses  are  therefore  less. 

As  the  burner  off-design  mass  flow  and  stagnation  values  have 
already  been  established,  the  off-design  conditions,  i.e., 
burner  inlet  Mach  number  and  recovery;  can  be  computed  once 
the  diffuser  recovery  and  variation  of  compressor  rotor  outlet 
Mach  number  M2;With  pressure  ratio  ffc,are  known.  Using  the 
dump  diffuser  model  and  the  variation  of  M2  with  it  shown 
in  Figure  14,  the  overall  diffuser  recovery  is  defined.  The 
results  are  shown  in  Figure  15  in  terms  of  burner  entrance 
Mach  number  M3  versus  percentage  of  power.  By  using  the  assumed 
combustion  model,  the  burner  recovery  is  evaluated  from 
Figure  5  at  this  value  of  M3  and  the  percentage  of  power.  The 
individual  diffuser  and  burner  recovery  along  with  the  com¬ 
bined  recovery  versus  percentage  of  power  are  shown  in  Figure  16. 


The  altitude  performance  was  determined  in  the  same  manner, 
with  the  assumption  that  M2  versus  ffc  is  invariant  with  alti¬ 
tude.  When  comparing  sea  level  to  altitude  performance  at 
the  same  value  of  ffc,  there  is  found  no  variation  in  M3  of  r) . 
However,  there  is  some  difference  in  M3  on  the  basis  of  per¬ 
centage  of  power  as  shown  in  Figure  15. 
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Figure  16.  Diffuser,  Burner,  and  Combined  Recovery  vs  Percentage  of  Power, Design 
Mach  Number  =  0.5. 


III.  BURNER  DESIGN 


This  section  describes  first  the  basic  concept  of  the  high¬ 
speed  burner  in  broad,  general  terms.  The  interrelations 
between  fuel  injection,  mixing,  burning,  pressure,  and  flow 
area  distributions  are  discussed  so  that  the  factors  leading 
to  the  chosen  configuration  may  be  understood.  After  this, 
the  three  fundamental  phenomena  associated  with  combustion 
chambers,  i.e.,  fuel  injection,  ignition,  and  combustion,  are 
discussed  in  detail. 

Only  liquid  JP-4  fuel  was  considered  for  the  following  reasons 
first,  to  eliminate  thermal  decomposition  and  deposition  pro¬ 
blems  accompanying  fuel  boilers;  and  second,  in  order  that 
the  high-speed  burner  concept  not  be  vitiated  by  the  con¬ 
sequence  of  specifying  only  a  fuel  in  the  vapor  phase.  As 
with  most  combustion  chambers,  the  present  design  is  compli¬ 
cated  by  the  necessity  of  providing  local  fuel-air  mixture 
ratios  that  are  near  stoichiometric.  Whereas,  in  conventional 
designs,  the  incoming  airflow  is  divided  to  form  a  primary 
combustor  zone  and  then  force-mixed  by  metering  holes;  in 
the  new  concept,  stoichiometric  zones  are  generated  spatially 
by  the  proper  injection  of  fuel.  Creating  the  zone  in  this 
manner  eliminates  the  enormous  total  pressure  loss  which 
accompanies  the  force-mixing  of  high-speed  airflows.  The 
method  and  analysis  used  to  produce  this  fuel  injection 
pattern,  along  with  other  problems  such  as  liquid  jet  break¬ 
up  and  droplet  heating  and  vaporization,  are  discussed 
under  the  heading  Fuel  Injection. 

The  combustion  and  mixing  processes,  besides  propagating  in 
the  spatially  varying  concentration  field,  must  also  be 
constrained  so  as  to  occur  at  constant  pressure.*  In  order 
to  satisfy  these  combustion  criteria,  the  GAST,  finite-rate 
combustion  with  turbulent  mixing  computer  program  was  used 
to  determine  the  requisite  burner  flow  area  distribution. 

A  number  of  computer  "experiments"  were  performed  in  the 


*  It  will  be  recalled  that  constant  pressure  combustion  was 
shown,  in  Section  II,  to  be  the  most  efficient  of  the  three 
types  of  burning  processes  examined. 


course  of  laying  out  a  satisfactory  burner  configuration. 
These  calculations  and  a  brief  description  of  the  computer 
program  are  given  under  the  heading  Combustion.  The  calcu¬ 
lations  showed  that  mixing  of  the  high-temperature-reacted 
products  with  the  cold  air  proceeded  very  slowly  and  was  the 
limiting  process.  This  phenomenon  necessitated  that  the 
characteristic  mixing  distance  be  the  smallest  dimension  - 
the  burner  flow  channel  height.  Thus  the  spatial  variation 
was  in  the  radial  direction,  with  the  flow  considered  to  be 
uniform  in  the  circumferential  coordinate. 


IGNITION 


Imposing  axial  symmetry  on  the  flow  necessitates  that  the 
burner  feature  circumferential  ignition.  Therefore,  the 
burner  design  included  a  combination  of  a  circumferential 
slot  flameholder  and  blunt-trailing-edge ,  spark-plug-driven 
strut-type  igniters  (called  primary  igniters)  .  Fuel  injected 
into  the  airflow  from  the  side  walls  of  each  primary  igniter 
is  induced  into  the  vortex  region  aft  of  the  blunt  trailing 
edge.  It  was  supposed  that  flame  established  initially 
in  the  wake  of  these  primary  igniters  would  subsequently 
spread  throughout  the  circumferential  slot.  In  this  way, 
the  main  fuel,  arriving  from  an  upstream  injection  station, 
is  ignited.  A  schematic  drawing  of  the  fuel  injection  and 
ignition  configuration  is  given  in  Figure  17.  It  can  be 
noted  in  the  figure  that  the  main  fuel  addition  is  at  station 
(F^ )  ,  the  igniter  fuel  is  at  station  (F2) ,  the  primary  igniter 
is  at  station  (P),  and  the  slot  flame-holder  is  at  station  (S)  . 


A  review  of  the  blunt-trailing-edge  flameholder  concept, 
given  in  Appendix  I,  shows  that  the  experimental  stability 
limits  of  this  igniter  can  be  correlated  in  terms  of  a  so- 
called  loading  parameter 


K(f/a)  = 


V 


,  .  .95  .85 

(p)  (D) 


(41) 


It  is  also  shown  in  Appendix  I  that  the  slot  flameholder 
stability  data  can  also  be  correlated  by  the  same  expression. 
Therefore,  both  the  slot  depth  and  the  trailing  edge  of  the 
primary  igniter  are  sized  according  to  Equation  (41). 
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FUEL  INJECTION 


The  problem  at  hand  is  to  inject  the  main  fuel  into  the 
burner  in  such  fashion  that  by  the  time  the  fuel  has  reached 
the  slot  flameholder,  it  is  vaporized  and  has  the  desired 
concentration  level  and  distribution.  These  criteria  form 
two  distinct  problems:  the  first  requires  knowledge  of  heating 
and  vaporization  rates  of  the  fuel  droplets  so  that  the  loca¬ 
tion  of  the  main  fuel  injector  relative  to  the  slot  can  be 
fixed;  the  second  requires  the  choice  of  the  proper  initial 
injection  pattern  and  the  knowledge  of  the  mixing  rates.  Only 
normal  injection  is  considered,  primarily  for  ease  of 
manufacture  and  installation. 

In  fixing  the  main  fuel  injector  location,  there  are  three 
distances  which  must  be  determined.  These  are 
(listed  in  order  from  the  injection  orifice): 

1.  liquid  jet  breakup  length,  x^ 

2.  length  required  to  heat  droplets  to  boiling 
temperature,  x 

H 

3.  length  to  vaporize  droplets, 

The  distance  between  injector  station  and  slot  is  then  taken 
as  the  sum  of  x^,  xH,  and  x^.  In  order  to  calculate  items  2 
and  3,  the  average  droplet  size  and  penetration  distance,  in 
addition  to  the  forces  on  the  droplets  and  th^ir  trajectories, 
are  determined  in  a  simple  manner. 

JET  BREAKUP 


When  a  liquid  jet  flows  out  of  an  orifice,  it  becomes  unstable 
due  to  effects  such  as  surface  tension,  and  breaks  up  into 
small  droplets.  If  the  liquid  is  injected  contracurrent 
or  normal  to  a  high-speed  airstream,  as  is  the  case  in  the 
present  burner,  high  shearing  stresses  are  set  up  on  the 
liquid  surface  and  cause  the  jet  to  shatter.  The  distance 
to  breakup  is  a  function  of  these  stresses, and  as  shown  by 
the  data  compilation  listed  in  Figure  2,  Reference  8,  the 
determining  parameters  are  the  Reynolds  number  and  Weber 
number. 
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In  the  present  case,  the  Reynolds  number  is 

pud 

N  =  L  —  =  4750 
Re  p 

L 

and  the  Weber  number  is  2 

p  V  d 

NWe  "  =  8600 


where 


(42) 


(43) 


UL  =  the  fuel  jet  velocity,  10  ft/sec 
V  =  air  velocity,  800  ft/sec 
PL  =  fuel  density,  48  lbm/ft3 
p  =  air  density,  0.32  lbm/ft3 

u 

A  *1 

PL  =  fuel  viscosity,  2.7  •  10~  lbm  ft  /sec 

CTL  =  fuel_air  surface  tension,  4.65  •  lo4  lbm/sec2 

03  2 

d  =  orifice  diameter,  - — —  ft 

w  J.  z 


so  that  from  the  aforementioned  data  compilation  chart,  this 
fuel  jet  corresponds  to  the  regime  called  immediate  jet 
breakup;  thus,  =  0. 

Penetration  Distance  and  Droplet  Size 


The  distance  that  the  fuel  penetrates  into  the  airstream, 
along  with  the  maximum  and  volumetric  mean  diameter  of  the 
resulting  fuel  droplet,  can  be  computed  with  the  above 
values  of  Reynolds  and  Weber  numbers.  Using  the  analysis  of 
Ingebo  in  Reference  9,  it  is  found  that  the  penetration 


distance  y^  equals 


y  =  d  (1.8) 

m  o 


N 


Re 


N, 


We 


0.7 


=  0.03  in. 


(44) 


the  maximum  drop  diameter  is 


D  =  d  (22.3) 
Mo 


VV 


v_ 

u. 


=  5  *  10 


-5 


ft; 


-0.29 


(45) 


and  the  volume-mean  drop  diameter  is 
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2.5  •  10"5  ft;  r  3Q  =  1.25  •  lo"5  ft 


Droplet  Trajectory 


The  following  simple  analysis  determines  the  distance 
necessary  to  accelerate  the  fuel  droplets.  It  is  assumed 
that  the  normal  velocity  components  can  be  ignored  (note 
the  negligible  penetration  distance) ,  that  all  droplets  are 
the  same  size  and  that  the  properties  of  the  oncoming  air¬ 
flow  are  constant. 

The  force  on  the  droplet  is  set  equal  to  the  drag  on  a  sphere 
so  that  from  Newton's  law, 


|pa  (V-U)2  CD  «  r20 


2  du 

r30  PL  U  dx 


(47) 


where  U  is  the  droplet  velocity,  is  the  drag  coefficient,  and 

x  is  the  axial  distance.  The  value  of  CD  is  a  function  of  the 
droplet  Reynolds  number,  where  this  Reynolds  number  is  defined 
as 


N 


0,  (V-U)  d30 


Re ,  D 


(48) 


During  the  acceleration  period;  Cp  varies  (c.f.,  Reference  10) 
from  0.9  at  U  equal  0  to  6.0  at  U  equal  0.9V.  By  integrating 
Equation  (47)  ir.  a  step-by-step  fashion,  assuming  Cp  constant 
at  its  average  value  during  each  x  step,  givcj  for  the  in¬ 


cremental  distance  x^. 


x .  = 

3 


1 


200  C 


log 


|  V~U  1  +  v  "  _J_  _  1  1 


Dj  l 


V-U  . 
1 


V-U 


V-U. 


(49) 


where 


C_ 


'Dj 

x 

3 

U 

i 


the  average  value  of  the  jraq 
coefficient  in  step  j 
the  incremental  distance 
droplet  velocity  at  the  beginning 
of  the  step 
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At  U=0.9V,  the  total  distance  traversed  by  the  droplets  is 
x  =  L  x.  =  0.18  in.  (50) 

j  3 

It  will  be  seen  later  that  this  distance  is  an  order  of 
magnitude  less  than  the  heat-up  or  vaporization  distance^o 
that  the  droplet  may  be  assumed  to  start  impulsively  at 
velocity  U  =  800  fps . 

HEAT-UP 

To  calculate  the  distance  required  to  vaporize  tne  fuel 
droplets,  the  heating  processes  are  imagined  to  be  separable 
into  two  distinct  entities:  a  heat-up  length  and  a  subse¬ 
quent  vaporization  length.  To  decouple  these  two  phenomena, 
it  is  assumed  that  the  thermal  conductivity  of  the  liquid  is 
large.  Therefore,  to  find  the  heat-up  time  and  associated 
droplet  travel  distance,  only  the  thermal  boundary  condi¬ 
tion  is  solved,  or 

4  3  AT 

3  »  r30  (pCp*L  If  '  hA'VV  (51> 

thus  ignoring  the  droplet's  spatial  temperature  variation, 
dT/dr.  This  concept  is  not  quite  true  during 
the  latter  part  of  the  heat-up  period,  when  the  droplet 
surface  has  reached  the  boiling  point  (with  attendant  vapori¬ 
zation)  ,  since  the  droplet  core  is  cooler  because  of  the 
finite  rate  of  heat  conduction.  The  error  introduced  by  this 
approximation  is  shown  to  be  quite  small  and  conservative. 

In  these  calculations  it  is  again  assumed  that  the  droplets 
are  spherical  and  are  all  the  same  size,  r3Q  =  1.25  •  10" ft. 
Also,  the  initial  fuel  temperature  is  T^^  =  520°R;  the  ambient 
air  temperature  is  T3  =  120O°R?  the  ambient  pressure  is  10  at¬ 
mospheres;  and  the  fuel's  boiling  point  is  taken  as  its 
initial  saturation  temperature  at  the  arbient  pressure,  or 
T3  =  900°R . 

Since  it  has  been  shown  that  the  relative  droplet  velocity 
(V-U)  is  small,  the  Stokes-flow  result  can  be  used  to  de¬ 
termine  the  heat  transfer  coefficient.  In  this  flow  regime, 
the  Nusselt  number  is 


I 
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2.0 


(52) 


Comparing  this  value  to  the  approximate  time,  previously 
determined  as  1.45  *  10~  seconds,  reveals  the  difference 
to  be  negligibly  small. 


Vaporization 


After  the  droplet  has  reached  boiling  temperature,  further 
convective  heating  vaporizes  the  liquid,  with  the  fluid  tempera¬ 
ture  remaining  constant.  The  rate  at  which  a  single  droplet 
evaporates  is 

-  I  pL 4frr2  ^7  =  h  4ffr2(T3-Ts)  (58) 

where  t  is  the  heat  of  evaporation  and  equals  150  BTU/lb. 
Substituting  the  Stokes-flow  relation  for  h  gives  the  follow¬ 
ing  vaporization  time  as  a  function  of  droplet  size: 


T 

V 


L 


2k, (T  -T  ) 

Ca  J  O 


J  d(r)2=  3 
r30 


(1- 


(59) 


The  time  for  complete  vaporization,  i.e., 
inspection: 

Ty  =  3  •  10  4  sec 


when  r  =  0,  is  by 


(60) 


and  the  vapo  ...zation  distance  is 


xy  =  V(TV)  =  2.4  in.  (61) 

The  total  distance  traveled  by  the  droplet  during  heat-up 
and  vaporization  isx=Xh+xv=3.85  inches ; therefore, 
the  main  fuel  injector  is  located  4.0  inches  from  the  slot 
f lameholder . 


Mixing 

The  second  part  of  the  fuel  injection  problem  is  to  obtain  a 
specified  fuel-air  distribution  in  the  vicinity  of  the  slot 
flameholder.  It  was  desired  that  the  fuel-air  ratio  be 
maintained  axisymmetric  with  a  stoichiometric  condition  ad¬ 
jacent  to  the  slot,  zero  f/a  at  the  opposite  wall;  in  addition 
the  distribution  should  have  a  step-like  variation.  The 
analysis  used  to  determine  the  mixing  of  the  injected  fuel 
with  the  oncoming  airflow  is  shown  in  Appendix  II. 
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By  assuming  that  the  fuel  injectors  act  like  po.nt  sources, 
the  fuel-ait  ratio  downstream  of  each  orifice  can  bo  expressed 


as 


r,  V  r  vr2 

f/a  -  -  — —  exp  -  - 

4TTM  D„X  H  .  4D  X 


(62) 


No  distinction  was  made  between  the  liquid  and  vaporized  fuel; 
however,  the  value  of  D^/V  was  taken  constant  at  0.0004  ft 
and  was  thus  more  representative  of  the  liquid.  The  axial 
location  (x)  was  4  inches  -  the  distance  computed  for  complete 
vaporization.  As  indicated  in  Appendix  II,  because  of  the 
linearity  of  the  basic  diffusion  equation,  the  superposition 
principle  allows  the  determination  of  the  fuel  concentration 
when  there  is  more  than  one  contributing  injector.  Finally, 
the  effect  of  the  bounding  annulus  walls  is  accounted  for  by 
the  assumption  of  mirror-like  fuel  reflection.  Thus,  the 
fuel  concentration  at  any  point  can  be  written  as 


«« 


f/a  =  E  (f/a)  . 

■i  ^ 


(63) 


and  expanding  to  include  the  effect  of  the  bounding  walls, 


f/a  =  E  (f/a) .  +  E  (f/a) .  (64) 

i  1  i  1 

Direct  Reflected 

The  results  of  the  fuel-air  mixing  computations  are  shown 
in  Figure  18,  which  is  a  representative  segment  of  the  slot  plane. 
The  listed  numerical  values  are  the  calculated  fuel-air 
ratios  and  the  phantom  lines  indicate  the  upstream  main  fuel 
injector  configuration.  (This  injector  design  is  discussed 
in  detail  in  the  next  section.) 

It  is  seen  from  the  figure  that  the  desired  concentration 
features  are  grossly  obtained  with  the  chosen  fuel  injector 
configuration.  This  distribution  could  have  been  improved 
if  another  row  of  injector  tubes  had  been  inserted;  however, 
the  added  complexity  and  drag  of  such  a  system  decided  against 
its  inclusion. 
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COMBUSTION 


The  GASL  finite-rate  reaction  and  turbulent-mixing  computer 
program  was  used  to  define  the  flow  processes  and  thus  the 
burner  area  variation  in  order  to  provide  the  desired  constant- 
pressure  burner  flow.  This  program  is  reviewed  in  some  de¬ 
tail  in  Appendix  III.  As  the  flow  channel  in  the  present 
burner  was  an  annular  passage  with  a  height  to  radius  ratio 
over  10,  the  2-D  analysis  was  utilized.  A  series  of  computer 
"experiments"  was  performed  in  order  to  lay  out  the  burner 
design.  Some  of  these  experiments  will  be  pointed  out  below. 
First  discussed,  however,  are  the  initial  conditions  and  the 
method  whereby  they  are  computed  -  a  necessary  prerequisite 
for  the  computer  program.  Next,  a  flow  model  is  postulated  for 
the  slot  flameholder  flow  followed  by  a  discussion  of 
computed  isotherms,  burner  flow  area  distribution,  and 
effect  of  pilot  temperature. 

Initial  Conditions 


Since  the  overall  fuel-air  ratio  is  much  less  than  stoichio¬ 
metric,  the  fuel  must  be  mixed  with  only  a  portion  of  the 
incoming  air.  This  zone  type  of  fuel  distribution  is  schema¬ 
tized  in  Figure  19,  and  three  initial  zones  are  seen:  a  pure 
air  region,  a  mixed  fuel-air  region, and  the  slot  flameholder 
region.  The  latter  two  zones  are  discussed  next. 

1.  Slot  Flameholder 

To  define  the  hot,  recirculating  flow  in  the 

slot, the  following  assumptions  are  made: 

a.  Pilot  gas  axial  velocity  is  taken  as  10 
percent  of  the  free-stream  value  (see 
discussion  below) . 

b.  There  are  no  circumferential  or  radial 
pressure  gradients  or  velocity  components. 

c.  Thermodynamic  properties  correspond  to 
those  resulting  from  JP-4  combustion  at 

f/a  =  0.0405  or  f/a  =  0.0676  (stoichiometric). 

d.  One-half  of  the  slot  depth  is  considered  in 
order  to  account  for  the  recirculating  flow. 
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Figure  19.  Definitions  Of  Coordinate  System. 
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Those  assun  pt  1  on coupled  with  the  sIol  dimension 
determine  the  pilot  muss  flow.  Previous  computa¬ 
tional  experience  at  GASL  has  shown  that  when  the 
level  of  piloi  to  froe-strcam  mass  flows  is  small 
(1/50)  changes  in  this  ratio  have  little  effect 
on  combustion  characteristics.  Thus  the  above 
procedure  can  be  used  to  define  the  initial 
pilot  flow  data  to  study  mainstream  flame  initia¬ 
tion  and  propagation. 


2 .  Fuel-Air  Region 


Since  the  fuel  is  injected  into  the  burner  in 
the  liquid  phase  and  at  a  lower  temperature  than 
the  airstream,  the  following  analysis  was  used 
to  determine  the  state  of  the  fuel-air  mixture. 

It  was  assumed  that  the  flow  was  one-dimensional 
and  of  constant  area  and  that  only  the  streamtube 
containing  the  fuel  had  contributed  to  its 
vapor ization  or  sensible  enthalpy  increase.  As 
the  fuel  was  injected  normal  to  the  airflow,  its 
momentum  flux  was  ignored.  Therefore,  from  con¬ 
servation  of  momentum. 


m  m 

pA"  +  — —  V  =  p  A"  +  —  V 
a  g  a  g  m 


(65) 


energy , 


Hm=  «  Ha  +  CHf+0=e  +  ^  (66) 

m  l+(f>  3  l+(f)  f  m  2 

a  a 


and  mass. 


•  •  •  •  f  -  £ 

m  =  m  +  nu  =  m  L  1  +  —  J 

m  a  f  a  a 


(67) 


along  with  the  perfect  gas  equation  of  state, 
R 


p  =  p  -  t  =p  +  p 

m  m  mw  m  a 
m 


(68) 


and  the  identity  concerning  the  flow  areas. 


A" 

_a 
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f /a 
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(69) 
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the  flow  parameters  »d  jno'nt  to  the  pilot  slot  >irc 
found.  In  the  above  eg '.at ions 


p  =  partial  pressure 
A"  flow  area  containing  fuel 
II  =  stagnation  enthalpy 
e  =  static  enthalpy 

(f/a)"  =  fuel  air  ratio  in  the  region  A" 

The  design  point  ummixod  conditions  which  are 
necessary  for  these  calculations  are  listed 
in  Table  III. 


TABLE  III.  UNMIXED 

AIR  AND 

FUEL  CONDITIONS 

Item 

JP-4  Fuel 

Velocity,  fps 

800 

10 

Stagnation  Temp.,  °R 

1200 

500 

Static  Pressure,  psia 

139 

139 

Molecular  Weight 

28.9 

128 

Besides  the  stoichiometric  fuel-air  layer, 
a  saturated  fuel  layer  and  a  fuel-air 
ratio  of  0.405  layer  were  also  evaluated.  The 
saturated  condition  was  the  extreme  case  where 
the  amount  of  fuel  just  necessary  to  saturate  a 
given  stream  tube  of  air  was  determiner*.  The 
smaller  f/a  case  showed  the  effect  of  lower  flame 
temperature  on  burner  design.  The  mixture  condi¬ 
tions  at  the  slot-pilot  for  these  fuel  layers 
are  given  in  Table  IV. 


# 
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Roflu  Its  of  Compul a  lions 

Tho  calculations  «»  -c  boqun  at  the  entrance  plane  of  the  slot. 
The  radial  distance  Y,  sec  Figure  10,  is  initially  at  x  =  0, 
measured  from  one- ha  1 f  the  depth  of  the  slot,  and  the 
distance  Yp  represents  the  'ntir-  flow  dimensions  and  is  thus 
a  function  of  axial  distance.  in  subsequent  paraqraphs,  the 
isotherms  generated  by  the  s co i ch iometr i c  fuel-air  layer  and 
the  saturated  fuel-air  layer  '/ill  be  described  in  detail, 
but  first  some  of  the  combustion  phenomena  arc  briefly  out¬ 
lined  to  make  clear  the  meaning  of  the  various  isotherm 
distributions. 

When  j  combustible  mixture  is  ignited,  two  basic  phenomen i 
occur.  First,  there  is  an  incubation  period  during  which 
no  sensible  heat  is  released  but  which  creates  the  free 
radicals  necessary  for  the  reaction.  This  is  called  the 
ignition-delay  time.  The  second  phenomenon  features  the 
large  heat  release  associated  with  flames  and  is  called 
the  reaction  time.  In  the  present  burner  configuration, 
there  is  a  third  phenomenon  due  to  the  fuel  and  air  distri¬ 
bution:  this  is  the  mixing  which  occurs  subsequent  to  the 
combustion. 

In  Figure  20,  the  isotherms  produced  by  the  saturated  fuel- 
air  case  are  portrayed.  A  very  long  ignition  delay  dis¬ 
tance  is  indicated,  approximately  3  inches,  which  is  due 
to  the  time  taken  for  the  diffusion  of  external  oxygen  and 
reaction  with  the  high-temperature  fuel.  After  ignition, 
the  flame  spreads  rather  rapidly  in  a  very  intense  zone  for 
an  axial  distance  of  about  b  inches.  At  this  point  the  slope 
of  the  isotherms  decreases  indicating  that  most  of  the  fuel 
is  burned  and  that  only  the  mixing  with  the  unburned  air  is 
continuing.  Here  the  flame  has  filled  only  one-half  the 
flow  passage  and  was  thus  deemed  to  be  unacceptable. 

The  stoichiometric  fuel-air  layer  generated  isotherms  are 
shown  next  ir.  Figure  2i.  The  ignition,  combustion,  and 
mixing  zones  are  seen  again;  but  now  at  an  axial  distance 
of  6  inches,  the  flame  has  filled  three-quarters  of  the 
flow  channel.  The  turbulent  mixing  process  completes  the 
flow  homogenization  by  x  =  1.1  ft.  This  latter  point  can 
be  seen  from  the  temperature  profiles  computed  at  various 
axial  stations  and  showrn  in  Figure  22;  by  x  =  1.07  ft.,  the 
pattern  factor  is  12  percent. 
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Temperature, 


Figure  22.  Temperature  Profiles  at  Various  Axial  stations. 
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To  examine  the  sensitivity  of  the  piloting  and  flame  speed 
due  to  pilot  temperature  and  initial  fuel-air  ratio,  the 
case  of  f/a  =  0.0405  (the  lean  flammibility  limit)  and  3000°F 
pilot  was  examined.  Results  are  shown  in  Figure  23  in  terms 
of  the  isotherm  spread.  It  should  be  noted  that  this  figure 
shows  only  the  initial  portion  of  the  combustion  zone.  These 
results  show  that  ignition  and  combustion  were  possible  at  a 
low  pilot  temperature  and  that  the  ignition  delay  time  is  of 
the  same  order  as  the  stoichiometric  layer.  Flow  area  varia¬ 
tion  is  shown  in  Figure  24,  which  also  portrays  the  difference 
in  the  two  types  of  fuel  layers.  It  may  be  seen  that  in  both 
cases  extensive  flow  area  variation  is  predicted  downstream 
of  the  flameholder  slot.  Since  the  stoichiometric  fuel-air 
layer  provided  the  smaller  burner  length,  this  configuration 
was  adopted  for  the  burner  design. 

Thus,  from  the  analysis  of  liquid  droplet  breakup  and 
organization,  fuel-air  mixing,  flame  stability,  and  the 
chemical  kinetics  and  turbulent  mixing  processes,  the  basic 
layout  of  the  burner  was  defined.  The  following  section 
discusses  the  burner  and  test  facility  design. 


I 


52 


i 


Pilot 


IV.  DESCRIPTION  OF  FACILITY 


Along  with  the  high-speed  burner,  this  section  describes 
the  test  facility,  the  associated  fuel  pumping  and  metering 
system,  and  the  burner  instrumentation.  Since  the  burner 
entrance  conditions  had  to  be  varied  over  a  rather  wide 
range  in  order  to  duplicate  the  altitude -power  spectrum 
shown  in  Section  II.  the  test  facility  featured  variable 
mass  flow,  stagnation  temperature,  and  pressure.  The 
technique  used  to  control  the  burner  entrance  flow  is 
described . 

A  schematic  of  the  test  facility  is  shown  in  Figure  25. 

TEST  SETUP 


The  burner  model  is  direct-connected  to  the  GASL  vertical 
pebble  bed  heater.  This  heater  consists  of  a  1500  psia 
maximum  pressure  vessel  containing  alumina  pebbles  heated 
by  "Globar"  electrical  resistance  heaters.  This  facility 
is  capable  of  providing  complete  simulation  of  the  actual 
engine  mass  flow,  temperature,  and  pressure. 

The  outlet  of  this  heater  is  a  sonic  orifice,  and  therefore 
the  airflow  sustains  subsequent  diffusion  to  subsonic  con¬ 
ditions  by  a  shock  system.  The  level  of  stagnation 
pressure  behind  the  shock  is  dependent  upon  the  shock 
location  so  that  pressure  levels  from  the  heater  stagna¬ 
tion  value  to  almost  ambient  can  be  obtained. 

From  the  pebble  bed  heater,  the  flow  enters  a  settling 
chamber,  which  consists  of  a  constant-area  duct  about 
14  inches  long  with  a  converging  section  at  the  downstream 
end  to  match  the  combustor  diameter.  Because  of  the  length 
of  the  settling  chamber  and  other  piping  upstream  of  the 
burner  and  the  relatively  narrow  burner  annulus  height,  the 
boundary  layer  was  calculated  to  be  too  large  and  would 
mask  the  high-speed  combustion  phenomena.  A  boundary  layer 
ram  scoop  is  located  at  the  leading  tip  of  the  converging 
section.  The  height  of  the  ram  scoop,  .070  inch,  is 
equal  to  the  calculated  local  boundary  layer  displacement 
thickness.  The  flow  rate  of  air  removed  from  the  facility 
is  continuously  measured  with  choked  venturis  so  that  the 
burner  inlet  airflow  is  precisely  known. 
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The  burner  is  formed  by  the  annulus  composed  of  a  cylindrical 
outer  shell  and  a  conical  centerbody.  The  outer  diameter 
of  the  burner  flow  section  is  5.05  inches  and  the  annular 
height  is  0.232  inch  at  the  igniter  plane.  The  burner 
measures  12  inches  from  the  slot  to  the  exit  and  the  cone 
half  angle  is  1°  29'.  The  resulting  flow  area  distribution 
is  different  from  the  one  prescribed  by  the  computer  analysis 
because  it  was  believed  to  be  more  prudent  to  test  first  a  simply 

shaped  body  having  the  design  J  ength  and  inlet  to  exit 
flow  area  ratio. 

The  burner  centerbody  which  is  shown  in  the  photograph, 

Figure  26,  is  supported  at  three  axial  locations  for 
structural  rigidity  and  concentricity.  These  locations 
are  at  the  model  leading  edge,  at  the  primary  igniters  and 
at  the  "Aero"-plug  exit.  The  photograph  also  shows  the 
boundary  layer  scoop  exhaust  tubes  as  well  as  the  burner 
attachment  flanges.  The  burner  section  was  designed  so  that 
it  could  easily  be  disassembled  for  modification  and  repair. 

Figure  2  7  shows  the  outer  shell  of  the  burner  mounted  to 
the  attachment  flange.  The  blunt-body-pilot  attachment 
bosses  and  the  rear  instrumentation  and  view  housing  are 
seen  in  the  figure. 

( 

A  drawing  of  the  ignition  region  is  shown  in  Figure  28, 
where  the  primary  igniters  and  the  circumferential  slot 
are  detailed.  To  be  noted  are  the  fuel  galleries  and  the  ^ 

spark  plug  located  in  the  blunt-trailing-edge  igniter. 

The  main  fuel  injector  configuration  that  was  chosen  is 

shown  in  Figure  29.  There  are  40  co-planar  injection 

orifices,  20  each  of  the  flush  and  the  off-wall  types, 

emanating  from  a  common  plenum.  Each  orifice  has  the  same 

outlet  flow  diameter:  0.032  inch.  The  off-wall  injector 

protrudes  approximately  0.10  inch  into  the  airflow.  * 

The  downstream  portion  of  the  test  burner's  outer  shell  is 
capable  of  rotating  360°.  This  arrangement  facilitates 
circumferential  surveying  of  the  annular  flow  field.  Two 
diametrically  opposed  Pyrex  windows  may  be  installed  in 
this  section.  Each  window  provides  a  2.25-inch  viewing  port 
and  is  used  for  observation  and  photography.  In  addition, 
the  ports  are  designed  to  be  used  as  mounts  for  instrumentation 
rakes.  Because  of  the  high  temperatures  prevailing  in  the 
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Burner  Centerbody. 
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Figure  27.  Viewing  Port  and  Instrument  Port  Housing. 


Spark  Plug 


(TYP) 


Figure  29.  Layout  Of  Main  Fuel  Injector  Configuration. 


combustor  flow  field,  the  rakes  are  recessed  in  a  cavity  and 
kept  out  of  the  flow  until  needed.  On  signal,  a  pneumatic 
actuator  inserts  and  withdraws  the  multiprobe  rakes. 

To  obtain  the  desired  pressure  level,  an  aerodynamic  plug 
valve  was  located  downstream  of  the  test  burner.  This 
plug  valve  consisted  of  a  geometric  constriction  and  varying 
amounts  of  cold  air  injected  to  create  a  sonic  condition 
at  the  plug.  Thus,  at  a  given  burner  mass  flow,  the  entrance 
pressure  level  is  varied  by  the  rate  at  which  cold  "plug" 
air  is  added.  The  plug  airflow  required  for  this  control 
system  is  shown  in  Figure  30  as  a  function  of  fuel-air  ratio 
for  burner  inlet  air  conditions  corresponding  to  the  100 
percent  and  60  percent  power  levels. 

The  design  of  a  combustion  test  model  included  a  heat  trans¬ 
fer  and  stress  analysis  and  the  evaluation  of  a  number  of 
thermal  protection  schemes.  These  are  detailed  in  Appendix  IV, 
where  it  is  shown  that  sufficient  heat  sink  cooling  is  avail¬ 
able  for  relatively  long  experimental  tests. 

I  NS TP UMENTAT I  ON 


The  instrumentation  rake  had  three  stagnation  pressure  and 
three  stagnation  temperature  probes  ,  fixed  side  by  side, 

1/8  inch  apart.  The  pitch  between  each  probe  was  1/8 
inch.  The  stagnation  temperature  probes  were  miniature 
unshielded/ceramo-type  platinum-platinum,  10  percent 
rhodium  thermocouples.  In  addition  to  the  rake  measurements, 
wall  static  pressure  and  wall  temperatures  were  measured. 

A  Scanivalve  pressure  distributor  was  used  to  monitor 
the  many  pressure  taps;  however,  continuously-reading  trans¬ 
ducers  were  used  to  monitor  the  more  important  pressures, 
such  as  the  burner  stagnation  pressure  and  fuel  pressure. 

Fuel  System 

A  schematic  diagram  of  the  fuel  system  is  shown  in  Figure  31. 
Three  independent  fuel  supplies  were  included  in  the  system. 
One  of  these  was  used  for  the  pilot  fuel  and  the  others 
supplied  the  main  fuel.  The  indicated  venturis  measured 
the  liquid  flow  rates,  and  since  they  were  a  cavitating 
type, they  also  uncoupled  the  flow  rates  from  any  downstream 
pressure  variations  thus  divorcing  the  fuel  system  from 
oscillations  in  the  combustion  chamber. 
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Figure  30.  Typical  Plug  Mass  Flow  Rates. 
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Figure  31. 


Fuel  Piping  Schematic  for  Army/GASL 
Combustor. 


The  fuel  system  was  calibrated  prior  to  testing  since  the 
cavitating  phenomenon  precludes  absolute  knowledge  of 
the  orifice  coefficients,  recoveries,  and  range  over 
which  the  venturis  provided  cavitating- type  performance. 

The  results  of  the  calibration  tests  are  shown  in  Figure  32 
for  two  different  venturis. 


Fuel  Pressure,  psia 
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Figure  32.  Calibration  Results  of  Cavitating  Venturis. 
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V. 


PRESENTATION  AND  DISCUSSION  OF  THE  EXPERIMENTAL  DATA 


It  was  desired  to  accomplish  the  following  tasks  in  the 
experimental  portion  of  the  program:  first,  to  demonstrate 
ignition  and  combustion  over  the  wide  range  of  operating  con¬ 
ditions  with  the  small,  annular,  high-speed  burner  con¬ 
figuration;  second,  to  obtain  predicted  performance  and  thus 
prove  the  viability  of  the  burner  concept;  and  finally  to 
improve  burner  performance,  e.g.,  by  eliminating  the  primary 
igniters.  Only  the  first  item  was  unequivocally 
demonstrated. 

No  difficulty  was  experienced  in  lighting-off  even  at  burner 
inlet  Mach  numbers  as  high  as  1.  Circumferentially,  uniform 
stagnation  pressure  and  temperature  distributions  could  not 
be  achieved,  with  the  result  that  approximately  25  percent  of 
the  total  volumetric  flow  was  devoid  of  combustion.  As 
these  unburnt  regions  are  located  centrally,  i.e.,  between 
primary  igniter  stations,  it  is  postulated  that  the  slot  pilot 
did  not  properly  initiate  combustion  throughout  the  entire 
circumferential  length. 

The  testing  procedure  called  for  mapping  the  entire 
operating  range  at  each  circumferential  station  before  pro¬ 
ceeding  to  the  next  station.  Since  the  test  conditions  of 
inlet  Mach  number,  total  temperature,  and  fuel -to -air  ratio 
could  not  be  kept  constant  from  run  to  run;  a  normalization 
procedure  was  adopted  to  portray  the  circumferential  distri¬ 
butions  of  the  significant  quantities.  With  respect  to  the 
temperature  measurements,  the  normalizing  procedure  was  to 
divide  the  measured  temperature  by  the  value  of  the  adiabatic 
flame  temperature,  Taf, where  Taf  was  determined  from  the 
measured  mass  flows  of  air  and  liquid  JP-4  and  the  GE  tables 
(Reference  3).  This  normalization  procedure  introduced 
other  problems,  which  are  discussed  under  the  heading 
Temperature  Distribution.  With  regard  to  pressure,  the 
recovery  is  plotted  which  effectively  normalizes  the  measured 
pressure  by  the  upstream  stagnation  value.  Figure  33  defines 
the  circumferential  stations  where  data  were  taken. 

The  nominal  flow  conditions  for  each  power  setting  are  list¬ 
ed  in  Table  V.  A  light- off  case  is  also  listed  which  demon¬ 
strated  the  capability  of  the  test  burner  to  initiate  com¬ 
bustion  under  simulated  startup  conditions. 


67 


TABLE  V.  INLET  CONDITIONS  FOR  ARMY/GASL  HIGH-SPEED 
COMBUSTION  TESTING 


Altitude 

(ft) 

Power  Setting 
(percent) 

m 

(lbm/sec) 

T3 

(  °R  ) 

P3 

(psia) 

f/a 

Sea  Level 

100 

5.0 

1190 

165 

.032 

25,000 

100 

2.0 

900 

61 

.024 

Sea  Level 

60 

4.37 

1120 

131 

.025 

25,000 

60 

1.65 

930 

45 

.02 

Sea  Level 

20 

2.96 

970 

82 

.016 

25,000 

20 

1.20 

795 

31 

.015 

Sea  Level 

Light-Off 

0.7 

520 

15.2 

.021 

The  data  are  portrayed  as  a  function  of  the  two  main  corre 
lating  parameters:  fuel-to-air  ratio  and  burner  inlet  Mach 
number.  These  two  parameters  are  determined  from  the 
experimentally  measured  values  as  follows; 

1.  The  fuel  flow  rate  is  fixed  by  the  venturi 
orifice  dimension  and  the  measured  nitrogen 
driving  pressure.  The  airflow  rate  is  the 
difference  between  the  measured  airflow  at 
the  exit  of  the  pebble  bed  heater  and  the 
measured  boundary  layer  bleed  flow. 

2.  The  burner  Mach  number  is  computed  from  the 
measured  inlet  air  stagnation  pressure  and  the 
average  of  seven  static  pressure  taps 
located  in  the  plane  of  the  primary  igniters 
just  upstream  of  the  blunt  trailing  edge. 

TEMPERATURE  DISTRIBUTION 


The  thermocouples  used  to  measure  the  total  temperature  had 
a  relatively  small  sensing  bead  diameter  compared  to  their 
lead-in  wire  dimensions.  Because  this  gives  rise  to 
large  conduction  errors,  an  analysis  was  undertaken  to 
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determine  the  necessary  correction  and  is  given  in  Appendix 

V.  Using  this  analysis,  a  temperature  correction  chart  is 
shown  in  Figure  51  with  which  the  true  gas  temperature  is 
determined  from  the  measured  junction  temperature. 

The  measured  temperatures,  corrected  according  to  Figure  51, 
are  shown  in  Figures  34  through  37  as  a  function  of 
circumferential  angle  0.  The  flow  conditions  associated 
with  each  graph  are  given  in  Tables  X  through  XIII  in  Appendix 

VI.  No  data  are  shown  for  the  20-percent  power  condition 
because  their  distribution  is  similar  to  Figure  36.  The 
burner  inlet  total  pressure  was  used  as  the  basic  flow 
property  to  define  the  relevant  power  setting.  No  attempt 
was  made  to  correlate  the  total  temperature  distributions 
according  to  burner  inlet  Mach  number  because  of  the  lack 
of  sufficient  data.  The  measured  circumferential  tempera¬ 
ture  distributions  reveal  the  following  pattern:  (1) 
downstream  of  the  igniters,  temperature  is  higher  than 
adiabatic  flame;  (2)  moving  away  in  the  0  direction,  the 
temperature  drops  below  adiabatic  and  remains  constant  at 
a  value  between  75  and  85  percent  of  T  f,  depending  upon 
power  setting;  and  finally,  (3)  midway  between  igniter 
stations,  combustion  terminates  very  rapidly. 

The  problems  of  associating  an  overall  fuel-air  ratio  with 
these  measurements  are  now  made  apparent.  The  addition  of 
primary  igniter  fuel  causes  the  fuel  distribution  to  be 
circumferentially  nonuniform;  and  because  of  the  low  mixing 
rates,  the  temperature  directly  downstream  of  the  igniter 
reflects  this  greater  f/a  and  thus  is  higher  than  the 
adiabatic  flame  temperature  computed  with  the  average  fuel- 
air  ratio.  If  the  actual  fuel-air  ratio  were  known  at 
each  point,  the  temperature  could  be  normalized  by  its  prop¬ 
er  value  so  that  the  combustion  efficiency  might  be  in¬ 
dicated.  Consider  the  following  proposed  fuel-air  distribu¬ 
tion  which  accounts  for  the  fuel  injected  from  the  primary 
igniters  (about  20  percent  of  the  total):  assume  that  one- 
half  this  fuel  is  diffused  out  by  0  =  45°  and  all  by  0  =  75°. 
Associating  values  of  Taf  to  these  equivalence  ratios  and 
normalizing  the  observed  temperatures  gives  the  dotted 
line  marked  "varying  f/a"  in  Figure  34.  If  this  circumfer¬ 
ential  fuel  distribution  is  as  postulated,  then  indicated 
combustion  efficiency  would  be  higher  than  that  implied  by 
the  normalization  based  on  a  global  T  value. 
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Stagnation  Temperature  Ratio 


Figure  34.  Circumferential  TemperaLure  Distribution, 
25,000  ftj  100  Percent  Power. 
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Figure  35.  Circumferential  Temperature  Distribution, 
25,000  ftj  60  Percent  Power. 
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Stagnation  Temperature  Ratio 


Figure  37.  Circumferential  Temperature  Distribution 
Sea  Level,  60  Percent  Power 


STATIC  PRESSURE  DISTRIBUTION 


A  typical  example  of  the  axial  distribution  of  static 
pressure  is  shown  in  Figure  38  where  the  pressure  is 
normalized  by  the  burner  inlet  stagnation  pressure.  The  line 
through  the  data  bars  represents  the  no-burning  distribution, 
while  the  data  points  indicate  the  pressures  during  combustion. 
These  data  are  from  two  different  runs,  both  with  identical 
burner  inlet  conditions.  It  is  seen  that  without  burning, 
the  static  pressure  increases  with  axial  distance  because 
of  the  increasing  burner  flow  area;  heat  addition  acceler¬ 
ates  the  flow  towards  Mach  1  so  that  these  two  effects  cause 
the  resulting  static  pressure  distribution  to  be  fairly  flat. 

In  addition,  there  is  practically  no  circumferential 
variation  in  the  static  pressure.  Another  axial  distri¬ 
bution  of  static  pressure  without  burning  is  shown  in  Figure 
39  for  the  case  of  an  inlet  Mach  number  of  0.58.  If  it  is 
assumed  that  all  the  total  pressure  loss  occurs  at  the  igni¬ 
ter  plane  and  the  pressure  recovery  is  then  one-dimensional, 
pressure  distributions  for  various  recoveries  may  be  computed; 
these  are  also  shown  in  Figure  39.  It  is  seen  that  a 
recovery  between  0.900  and  0.915  just  brackets  the  data 
and  serves  to  indicate  the  average  loss  due  to  the  dump 
diffusion  process. 

total  pressure 


It  was  found  that  the  pressure  recovery  was  a  function  of 
the  burner  inlet  Mach  number,  fuel-to-air  ratio  and 
circumferential  location.  The  pressure  recoveries  at  three 
circumferential  stations  (0°,  30  and  90°)  at  various  fuel- 
air  ratios  are  shown  in  Figures  40,  41, and  42  respectively. 
It  may  be  seen  that  there  is  little  or  no  change  in  the 
results  due  to  variation  in  pressure  level.  The  circum¬ 
ferential  pressure  recovery  distributions  are  obtained 
from  the  foregoing  charts  and  from  other  data  at  45°,  75° 
and  225°.  At  the  design  Mach  number,  this  distribution 
is  shown  in  Figure  43.  Because  of  geometrical  symmetry 
about  the  primary  igniter  centerline,  it  was  assumed  that 
the  0°  data  point  could  be  reflected  to  0  =  60°.  Note 
that  reflecting  the  0  =  2  2  5°  data  point  to  0  =  75°  shows 
fairly  good  agreement  with  the  actual  measurement  taken. 

The  effect  of  Mach  number  on  the  circumferential  pressure 
distribution  is  shown  in  Figure  44  at  a  constant  fuel-air 
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^^ithout  Burning 


Figure  38.  Axial  Static  Pressure  Distribution  With 

and  Without  Burning,  m3=0.48,  P  =128  psia 


Static  Pressure,  p/P 


Figure  39.  Pxial  Distribution  of  Static  Pressure, 


M  =  0.58,  No  Burning- 
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Figure  41.  Burner  Pressure  Recovery  vs  Mach  Number 
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rcumferential  Angle,  ©^Degrees 


Effect  of  Burner  Mach  Number  on  Pressure  Recovery,  f/a  =  0.018  . 


ratio  of  0.018.  It  appears  that  the  shape  of  the  distri¬ 
bution  is  unchanged  by  the  Mach  number,  and  only  the  level 
of  the  recovery  is  reduced  by  increasing  the  Mach  number. 

COMPARISON  OF  EXPERIMENT  AND  THEORY 

As  the  experimentally  observed  burner  static  pressure  level 
was  fairly  constant  during  combustion,  a  comparison  of 
burner  recovery  can  be  made  between  constant  pressure  theory 
and  experiment.  Figure  45  shows  this  comparison  where  the 
experimental  data  were  all  taken  at  the  0  =  30°  station. 

It  is  seen  that  theory  fairly  well  predicts  the  observed 
recovery  especially  at  the  lower  values  of  fuel-air  ratio. 
The  agreement  is  particularly  interesting  since  the 
theoretical  recovery  values  do  not  include  any  pressure  loss 
due  to  flow  blockage  caused  by  the  primary  igniters  (Note 
the  value  of  at  f/a  =  CL)  A  possible  explanation  for 
this  mitigating  effect  is  that,  due  to  burning  in  the 
base  region  behind  primary  igniters,  there  is  a  reduction  in 
the  large  dissipative,  vortical  flow  usually  associated 
with  the  process  of  dumped  diffusion. 
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VI .  CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS 

From  an  analytical  and  experimental  program  studying  the 
high-speed  burner  concept,  the  following  items  are  deemed 
significant : 

1.  Consideration  of  the  combined  dump-diffuser  and 
constant-pressure  burner  efficiences  indicate 

an  optimum  burner  inlet  Mach  number  equal  to  0.5. 

2.  The  GASL  finite-rate  combustion  and  turbulent 
mixing  programs  predicted  the  completr  burning 
and  mixing  within  the  relatively  short  length 
of  1.1  ft  for  the  specified  geometry  and  fuel- 
air  layer  model. 

3.  Ignition  and  combustion  of  liquid  JP-4  fuel 
were  experimentally  obtained  over  the  range  of 
operating  conditions. 

4.  Combustion  was  manifest  in  the  vicinity  of 
the  blunt-trailing-edge  primary  igniters,  but 
not  in  the  region  between  igniter  stations.  It 
appeared  that  the  slot  flameholder  was  in¬ 
effective,  causing  some  25  percent  of  the  total 
fuel  flow  to  remain  unburned. 

5.  Comparison  of  experirrent  and  recovery  downstream 
of  the  primary  igniters  with  predicted  constant- 
pressure  combustion  recovery  shows  fairly  good 
agreement. 

RECOMMENDATIONS 

It  is  believed  that  the  reported  experimental  results  prove  the 
feasibility  of  the  high-speed  burner  concept .  However,  it  is 
clear  that  the  present  knowledge  of  the  working  of  the  slot 
flameholder  is  inadequate.  Therefore,  it  is  recommended  that 
future  effort  be  directed  toward  understanding  the  following 
flameholder  phenomena: 

1.  blow-off  stability  limits 

2.  mechanism  of  flama  initiation 

3.  effects  of  slot  geometry 
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APPENDIX  I 

FLAME  STABILIZATION  AND  IGNITION 


In  developing  high-output  combustors,  it  is  necessary  to 
stabilize  the  flame  in  a  mixture  of  fuel  and  high-velocity 
air.  This  stabilization  can  be  accomplished  by  placing  a 
bluff  object  in  the  gas  stream, and  if  conditions  are  chosen 
properly,  a  flame  will  be  stabilized  in  the  recirculation 
zone  formed  behind  the  bluff  body. 

A  number  of  investigators  have  determined  empirically  the 
range  of  air-fuel  ratios  for  which  bluff  bodies  of  various 
shapes  and  sizes  will  stabilize  a  flame  as  a  function  of 
mixture,  velocity,  pressure,  and  temperature.  These  experi¬ 
ments  were  somewhat  idealized  by  using  a  vaporized  and 
premixed  fuel-air  mixture.  This  is  important  because,  as  shown 
in  Reference  13,  flame  stabilization  is  affected  by  the 
amount  of  liquid  left  in  the  vaporizing  stream  and  drop 
size. 

The  nature  of  the  flow  in  the  recirculating  zone  of  a  flame 
stabilizer  was  illustrated  by  the  work  of  Nicholson  and 
Field  (Reference  14)  .  Small  particles  of  sodium  acetate 
were  placed  in  the  airstream.  When  one  of  these  particles 
entered  the  flame  zone,  the  resulting  sodium  vapor  emitted 
light,  which  was  recorded  by  high-speed  motion  pictures 
(1/5000  second).  These  motion  pictures  have  shown  very 
rapid  mixing  within  the  recirculating  region,  and  it  is 
generally  felt  that  stabilization  is  accomplished  by  ignition 
of  fresh  mixture  by  mixing  with  the  hot  gases  in  the  re¬ 
circulating  region.  The  hot  combustion  products  leaving 
this  zone  cause  ignition  of  the  main  body  of  fuel  and  air. 

A  quantitative  understanding  of  combustion  in  this 
important  eddy  region  would  require  much  more  detailed  in¬ 
formation  on  the  heat  balance,  local  combustion  efficiencies, 
flow  patterns,  etc.,  than  is  now  available.  Nevertheless, 
several  investigators,  by  making  assumptions  as  to  the  mechanism 
of  flame  blowoff,  have  attempted  to  establish  the  groups 
of  variables  expected  to  control  flame  stabilization 
(References  14,  15,  and  16) .  These  analyses  generally  arrive 
at  the  conclusion  that  the  mixture  ratio  at  blowoff  is  a 
function  of  v /paD^Tc,  where  v  is  the  velocity  past  the  bluff 
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body,  P  is  the  absolute  pressure,  D  is  a  characteristic 
dimension  of  the  stabilizer,  and  T  is  the  mixture  absolute 
temperature.  Figure  46  shows  a  typical  correlation  as 
published  by  DeZubay  (Reference  16),  showing  the  effect  of 
body  diameter,  pressure  and  mixture  velocity.  The  gas 
velocity  varied  from  40  to  550  ft/sec,  the  disc  diameter 
varied  from  1/4  in.  to  1  in.,  and  P,  the  absolute  pressure, 
varied  from  3  to  15  psia.  The  gas  was  a  propane-air 
mixture  initially  at  room  temperature.  The  fuel-air  ratio 
at  blowoff  is  plotted  as  a  function  of  the  dimensional 
group  V/p- #  where  D  is  the  diameter  of  the  circular 
flat  disc  used  to  stabilize  the  flame.  It  is  seen  that 
larger  bodies,  higher  pressures,  and  lower  velocities 
give  wider  stability  limits  and  that  a  maximum  value  of 
the  above  group  for  stable  operation  exists  in  the  region 
of  stoichiometric. 

In  the  work  of  Reference  17,  the  effect  of  pressure  was 
studied  and  was  also  found  to  have  an  approximately  inverse 
dependence  in  the  stability  criterion.  This  pressure 
effect  is  of  great  importance  for  the  application  of  interest 
here.  Flammability  limits  and  laminar  burning  velocity 
have  been  found  to  be  less  significant,  while  spark  ignition 
energy  may  have  a  larger  effect.  Similar  correlations 
of  the  effect  of  diameter  and  velocity  have  been  published 
by  other  investigators  (References  15,  18,  19). 

There  is  considerable  variation  in  the  exponent  of  D  (.45 
to  1.00)  and  it  is  believed  that  these  differences  are  too 
large  to  be  attributed  to  experimental  errors.  Table  VI 
summarizes  the  work  done  by  several  investigators.  The 
first  two  arrived  at  an  exponent  of  1.00  and  .86  using 
cones  or  flat  discs.  The  next  two,  using  rods  mounted  with 
their  axes  90°  to  the  direction  of  flow,  found  that  an 
exponent  of  .5  best  fitted  their  data.  It  is  apparent  from 
these  data  that  the  correlation  of  blowoff  as  a  function 
of  Da  is  not  completely  satisfactory.  It  would  appear  that 
a  more  accurate  study  of  flameholders  along  with  more 
knowledge  of  the  aerodynamics  as  a  function  of  velocity, 
shape,  and  heat  release  is  needed. 

The  effect  of  temperature  has  been  investigated  in  References 
18  and  19.  An  exponent  of  1.2  on  the  absolute  temperature 
T  best  correlated  the  data. 
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is  the  difference  between  outer  and  inner  radii;  for  rods  or 
is  the  width  measured  perpendicular  to  the  axis  of  gutter. 


Slot  type^^flal^holder^^1^6111191  points  for  th® 
of  Reference  20  it  i «  ich  was  computed  from  the  data 

for  *  lean'fuel-a^"868  ^^-^2  at^sT 

line.  RichUmixtures^aneLovi  theMnflnd0"^^  DeZubaV 

L?Lrg:?erateS  a  —  f lameholder^tha^a^luf f at 


Fuel  Air  Ratio  ,  f/a 


Figure  46.  Flameholder  Stability  Limits. 


93 


APPENDIX  II 

REVIEW  OF  MIXING  ANALYSI S 


The  rate  at  which  fuel  is  transported  radially  by  turbu¬ 
lent  eddy  diffusion  depends  on  the  time-averaged  fuel  con¬ 
centration  gradients  and  on  the  eddy  diffusivity  of  mass, 
characterized  by  the  turbulent  diffusion  coefficient 
D,p  (ft^/ sec)  .  For  axially  symmetric  flows,  the  turbulent 
diffusion  equation  may  be  written  as: 


where 
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y  =  local  mass  fraction  of  fuel  =  f/f+a  sa  f/a 
p  =  total  density  of  mixture 
v  =  radial  velocity  component 
u  =  axial  velocity  component 
r  =  radial  distance 
x  =  axial  distance 


In  writing  this  equation,  it  was  assumed  that  diffusion  in 
the  axial  direction  is  small  compared  to  convective  bulk 
flow.  Equation  (70)  may  be  rewritten  in  the  form 
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It  is  convenient  (although  not  accurate)  to  assume  that  Dt 
is  only  a  function  of  x  or  even  constant  in  a  given  airstream. 
In  addition,  if  compressibility  effects  are  neglected  and 
the  radial  velocity  is  very  small,  then  Equation  (71) 
becomes 
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The  solution  to  Equation  (72)  may  be  obtained  by  using 
standard  Fourier  transform  techniques  and  is  a  function  of 
the  boundary  and  initial  conditions.  The  solutions  for 
three  common  injector  configurations  are  shown  in  Table  VII 
taken  from  Reference  12. 


TABLE  VII.  SOLUTIONS  TO  THE  DIFFUSION  EQUATION 


Type  of  Fuel  Injector 


Solution 


Point  Source 


Disk  Source 
(of  Radius  R) 

Ring  Source 
(Ring  Radius  R) 
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where 


M  =  mass  flux  of  air  lb/sec,  ft 
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I  is  the  modified  Bessel  function  of  the  first  kind,  zero 

O 

order.  Both  <p  and  $  are  plotted  in  Reference  12  over  a 
range  of  flow  conditions. 
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Since  the  governing  diffusion  equation  is  linear,  the 
fuel-air  distributions  resulting  from  multi-source  in¬ 
jections  can  be  expressed  by  suranaries  of  the  individual 
source  equations.  This  means  that  the  fuel  in  a  volume  of 
air  is  a  summation  of  the  fuel  contributed  by  each  in¬ 
jection  source  acting  as  though  it  were  the  only  source 
present.  Therefore,  we  may  write  that 

7=2  (f/a)  (79) 

a  L  i 

The  solutions  described  to  this  point  are  subject  to  the 
limitation  of  infinite  duct  size.  Presence  of  a  wall 
implies  a  boundary  condition  of  the  type 
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and  consequently  the  general  solution  is  no  longer  valid. 


For  a  single  point  source  or  disk  source  located  axially  in 
a  finite  duct,  these  are  available  solutions  for  the  fuel 
distribution.  Thus  for  a  point  source  injecting  into  a 
finite  duct.  Reference  12  gives  the  following  expression: 
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where 


(f/a)  *  after  mixing  fuel-air  ratio 

o 

j  -  Bessel  function  ,  zero  order 
o 

R  =  duct  radius 

«  roots  of  J,  first-order  Bessel  function 
n  1 


Similar  expressions  can  be  derived  for  non-point  and  non- 
axial  fuel  injectors,  but  the  resulting  equations  become 
very  complex.  Therefore,  for  the  present  combustion  design, 
it  was  decided  to  use  a  graphical  technique  developed  by 
Longwell  and  Weiss  and  shown  to  be  a  good  approximation 
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in  Reference  12.  This  graphical  technique  in  based  on  the 
assumption  that  fuel  diffuses  from  the  wall  as  a  mirror 
reflection  of  the  fuel  which  would  have  diffused  beyond  the 
wall  if  the  wall  had  not  been  there. 


Turbulent  Transport  Properties 

Predicting  the  fuel  distribution  accurately  depends  upor 
the  utilization  of  the  correct  values  for  the  turbulent 
transport  coefficient.  The  terms  Dj  and  v  which  appear  in 
the  diffusion  equation  as  the  ratio  Dj/v  may  be  reckoned  as 
a  diffusion  parameter  which  has  the  dimension  of  length. 
Since  the  diffusivity  D-p  is  proportional  to  the  turbulent 
intensity  of  the  airstream,  one  would  expect  D-j>  to  vary 
airectly  with  the  flow  velocity  so  that  the  ratio  Dt/v 
should  be  constant. 

Longwell  and  Weiss  experimentally  determined  and  showed  in 
Rcierence  12  the  variation  of  D^/v  as  a  function  of  v  for 
a  very  volatile  liquid-naptha  and  for  the  diesel  fuel. 
Because  of  the  presence  of  the  liquid  droplets,  it  was 
found  that  the  value  of  D^/v  dropped  off  with  increasing  v . 
However,  the  value  for  naptha  approached  a  constant 
(Dt/v  «  .001)  beyond  v  =  400  ft/sec,  since  it  was  almost 
completely  evaporated.  Extrapolat  ing  the  diesel  fuel 
(kerosene)  value  to  the  air  velocity  of  the  present 
burner,  800  ft/sec,  gives  a  value  of  D<p/v  =  .0004  ft. 
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APPENDIX  III 

REVIEW  OF  FINITE-RATE  CHEMICAL  KINETICS 


For  an  analytic  description  of  high-speed  combustion 
problems,  it  is  generally  not  adequate  to  assume  equilibrium 
burning,  since  the  time  for  the  combustion  process  is  too 
short.  Hence,  it  is  necessary  to  formulate  and  employ  a 
finite-rate  chemical  kinetic  combustion  mechanism  for  the 
chemical  system  of  interest. 

Work  at  GASL,  as  described  in  Reference  2  ^resulted  in  the 
formulation  of  a  finite-rate  chemical  kinetic  model  whose 
numerical  solution  technique  was,  without  significant  loss 
in  accuracy,  at  least  two  orders  of  magnitude  faster  than 
standard  techniques  such  as  Runge-Kutta.  The  numerical 
solution  technique  of  Reference  21  was  improved  upon  in 
Reference  22,  which  (besides  decreasing  numerical  solution 
time  even  more)  made  it  practical  to  perform  finite-rate 
chemical  calculations  for  more  complex  systems,  such  as  the 
burning  of  a  hvdrocarbon-air  mixture. 

One  very  important  family  of  hydrocarbons  is  the  paraffins 
(CnH2n+2) •  Work  was  done  at  GASL,  as  reported  in  Reference 
23,  to  formulate  a  detailed  chemical  kinetic  model  for  the 
oxidation  of  the  paraffin  propane.  This  was  done  because 
the  oxidation  processes  inherent  in  the  propane  system  are 
quite  similar  to  heavier  paraffins.  As  borne  out  by 
experimental  results,  the  ignition  delay  and  reaction  times 
of  paraffins  even  heavier  than  kerosene  (C9H20)  are  within 
an  order  of  magnitude  of  those  for  propane.  The 
analysis  employed  31  chemical  species  participating  in  69 
elementary  reactions.  This  detailed  model  was  employed 
in  one-dimensional  stream  tube  type  calculations  over 
a  wide  range  of  pressures,  temperatures,  and  equivalence 
ratios.  Results  were  compared  with  data  from  both  shock 
tube  and  steady  flow  facilities  (References  24  and  25) . 

In  general,  the  computed  ignition  delay  times  were  well 
within  the  order  of  magnitude  of  the  experimental  data. 

However,  because  this  reaction  mechanism  for  the  hydrocarbon 
fuel  system  is  extremely  complex,  it  was  found  to  be  too 
expensive  for  use  in  flow-field  calculations.  In  addition. 
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the  lack  of  well-defined  rate  constants  for  many  of  the 
intermediate  reactions  justifies  seeking  smaller,  yet 
representative,  chemical  systems  for  the  combustion  of 
kerosene.  Therefore,  the  detailed  kinetic  model  was  used 
to  formulate  a  relatively  simple  "quasi-global "  kinetic 
model  that  would  be  appropriate  for  both  fuel  lean  and 
fuel  rich  paraffin-air  mixtures.  It  was  desired  to  pre¬ 
dict  both  ignition  delay  times  and  total  reaction  times  to 
within  an  order  of  magnitude.  The  use  of  a  "quasi-global” 
model  with  less  than  a  dozen  chemical  species  makes  it 
economically  feasible  to  perform  multidimensional  flow- 
field  calculations. 

The  species  employed  in  the  quasi-global  model  are  the 
paraffin  CnH2n+2  and  the  equilibrium  combustion  products: 
H2»  02*  h20,  c°2'  c0'  h*  0*  0H-  In  addition^N2  is  present 
as  a  diluent, and  solid  carbon  (graphite) may  also  be  con¬ 
sidered  to  be  present  in  rich  mixtures.  The  reaction 
mechanism  is  comprised  of  the  nine  intermediate  reactions 
shown  in  Table  VIII. 

The  production  rate  of  species  i,  in  moles/cc/sec,  is 
Qj”en  by 


18 
=  Z 
r*l 


{V"  -v\  )k 

i  ,  e  i ,  r  r 


9 

D  (C. )»  •  . 

j.l  =>  3 


,  r 


(82) 


for  the  reaction 


Z 

i 


i  l 


R. 

l 


(83) 


til 

where  Cj  is  the  molar  concentration  of  the  j  species, 
is  tne  stoichiometric  coefficient  i  associa  ted  with  the 
species  and  kr  is  either  the  forward  or  the  reverse  reaction 
rate  constant.  The  above  system  has  been  solved  by  the 
linearization  technique  described  in  Reference  21.  This 
technique  involves  the  expansion  of  the  product  term  in 
Equation  (82)  in  Taylor’s  series  and  truncating  to  the 
lowest  order  term. 
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TABLE  VITI.  REACTION  MECHANISM 


The  reaction  rate  constants  are  expressed  in  the  standard 
Arrhenius  form  given  by 

a  ^r 

kr(t)  =  At  r  exp  (  -  —  )  (84) 

where  A,  and  are  the  constants  given  in  Table  VIII. 
The  species  enthalpy-temperature  data  have  been  curve 
fitted  from  the  data  in  the  JANAF  Tables. 

The  quasi-global  reaction  was  written  as 


CnH2n+2  +  2  °2  nC0  +  <n+1)H2  (85) 


This  was  combined  with  the  above  nine  elementary  reactions 
for  the  interaction  of  the  more  stable  combustion  products. 
The  quasi-global  fuel  oxidation  reaction  rate  is  of  the 
f  orm 


dc 


C  H 
n  2n+2 

dr 


=  k  f(t) 


a- 

P 


C 


fi- 

CnH2n+2 


(86) 


The  activation  energy  E  was  based  on  experimental  data  and 
was  determined  to  be  13,740  cal/mole.  fi'  was  set  equal 
to  1/2  and  y'  equal  to  1.  Comparisons  with  the  detailed 
kinetic  model  for  different  temperature  pressures,  and 
equivalence  ratios  were  used  to  determine  the  values  of 
k,  f(t),  and  a'  .  The  resulting  fuel  burning  rate  is 


dC 


C  H  _ 
n  n+2 

dr 


=  2 


10c 


t(°K) 

1111 


-•5  I 


p (atm) 


.3 


(87) 


'C  H-> 

n  2n+2  0. 


13,740 

Rt 


The  applicable  temperature  range  of  the  above  quasi-global 
model  and  the  nine  elementary  reactions  used  in  conjunction 
with  it  is  in  the  region  of  800°  to  3000°K.  Comparison  of 
the  quasi-global  model  with  the  detailed  kinetics  model  of 
Reference  23,  and  with  experiments  for  ignition  delay  and 
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reaction  times,  results  in  agreement  to  within  an  order 
of  magnitude. 

Previous  work  at  GASL  (References  2  7,  28  and  29)  resulted 
in  the  formulation  of  a  numerical  finite  difference  technique 
for  the  solution  of  multi-dimensional  turbulent  mixing 
problems  employing  the  boundary  layer  equations  for  the 
general  case  of  arbitrary  Prandtl  and  Lewis  numbers  and 
axial  pressure  gradient.  Semi-empirical  eddy  transport 
coefficients  have  been  used,  based  on  experimental  results, 
such  as  the  work  presented  in  Reference  28. 

The  mixing  analyses  coupled  to  the  finite  rate  chemical 
system  have  been  used  to  analyze  axisymmetric  and  plane 
two-dimensional  hypersonic  wakes  (Reference  27),  free 
jets  (Reference  28),  and  ducted  configurations  (Reference 
29). 

The  combination  of  the  quasi-global  chemical  model  coupled 
to  the  finite-difference  mixing  analysis  is  the  analytical 
tool  used  for  designing  the  present  high-speed  combustion 
chamber. 
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APPENDIX  IV 

HEAT  TRANSFER  ANALYSIS 


A  heating  analysis  is  made  of  the  subject  combustor,  and 
a  number  of  thermal  protection  schemes  are  evaluated.  The 
heat  flux  is  calculated  by 


=  -O23 

W=  -.2 


i4 


.2 


2/3 
Pr  J 


N. 


<pv)^8  <T4-Tw)  (88) 


The  results  are  plotted  in  Figure  47  which  shows  the  effect 
of  combustion  gas  temperature  on  heat  flux  with  wall  temp¬ 
erature  as  parameters.  Note  that  with  combustion  tempera¬ 
tures  of  3000°R,  the  heat  flux  varies  from  300  BTU/sec-ft2 
for  a  cold  wall  (60°F)  down  to  130  BTU/sec-ft2  for  a  wall 
at  2000°R. 


Air,  fuel,  and  heat- sink  cooling  methods  were  evaluated, 
and  comparisons  were  made  to  determine  the  technique  that 
is  most  advantageous  for  this  application. 

Air  Cooling 

Air-cooling  requirements  were  calculated  using  a  heat 
transfer  equation  similar  to  that  presented  above,  and 
the  results  are  presented  in  Figure  48.  Here,  the  effect 
of  heat  flux  on  coolant-side  wall  temperature  is  shown 
with  coolant  air  velocity  as  a  parameter.  An  anticipated 
requirement  might  be  to  use  a  coolant  air  velocity  of  80 
ft/sec  at  150  psia  to  maintain  the  coolant-side  wall 
temperature  at  1600°F  while  subjected  to  a  heat  flux  of 
130  BTU/sec-ft  .  The  coolant  airflow  rate  would  then  be 
2.6  lbm/sec  for  a  1/8  inch  annulus  circumscribing  a 
6 -inch  diameter.  Air  cooling  for  the  centerbody  would 
have  similar  flow  requirements.  Maximum  bulk  heating  of 
the  coolant  air  would  be  of  the  order  of  200°  to  300°F. 


Fuel  Cooling 

Fuel-cooling  requirements  were  calculated  using 
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Figure  48. 


Effect  of  Heat  Flux  on  Coolant- 
Side  Wall  Temperature. 


(89) 


q 


conv 


ft  k  Pr 

,027  (pV)g  (Tw-Tg)(  — 8-  )g( 


1.4 
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and  liquid  JP-4  as  the  coolant.  Figure  49  shows  the  effect 
of  heat  flux  on  coolant  fuel-side  temperatures  with  fuel 
velocity  as  a  parameter.  A  total  fuel  flow  of  approximately 
0.2  lb/sec  corresponds  to  airflow  of  5  lb/sec  and  has 
only  10  percent  of  the  heat  capacity  required  to  cool  the 
entire  combustor.  Fuel  cooling,  therefore,  can  be  used 
for  cooling  specific  local  areas  but  not  the  entire  center- 
body  or  outer  liner  surface  areas. 

Heat-Sink  Cooling 

Heat-sink  cooling  characteristics  were  determined  by  using 
Schneider  charts  representing  analytical  solutions  for 
transient,  one-dimensional  heat  conduction  in  single— layer 
plates  with  constant  ambient  temperature  and  with  finite 
internal  and  surface  conductance.  Figure  50  3hows  the 
temperature  response  of  the  combustor  walls  with  complete 
combustion  as  a  function  of  time  with  wall  thickness  as  a 
parameter.  Figure  51  shows  the  effect  of  heat  flux  on 
wall  temperature  for  a  test  duration  of  60  seconds  with  wall 
thickness  as  a  parameter. 

It  can  be  observed  f’-om  these  figures  that  heat-sink 
cooling  subjects  the  vails  to  temperature  response  throughout 
the  entire  test  run  ard  to  very  high  metal  temperature 
exposure  during  thf.  latter  part  of  a  test  run;  however,  this 
technique  is  preferred  on  the  basis  of:  (1)  achieving 
minimum  required  wall  temperature  within  a  few  seconds, (2) 
simplicity,  and  (3)  lowest  cost. 

Structural  Considerations 


Thermal  stress  resulting  from  temperature  differentials 
across  the  combustor  liner  walls  may  be  determined 
approximately  by 


= 

t  2(1-0) 


(90) 
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r  ■  60  seconds. 


where 


E  =  the  modulus  of  elasticity 
6  =  coefficient  of  thermal  expansion 
£T  =  the  temperature  differential 
p  =  Poisson's  ratio 

Thin-walled  tube-type  wall  structures  minimize  thermal 

stresses;  however,  a  practical  shell-type  wall  thickness 

is  1/4  inch.  An  air-cooled  stainless  steel  wall  of  this 

o 

thickness,  when  subjected  to  the  2500  F  ambient  temperature 
gas  with  corresponding  fuel  flux  of  130  BTU/sec-ft2  and 
wall  temperature  differential  of  800°F,  sustains  thermal 
stress  amplitudes  of  140,000  psi. 


For  a  stainless  steel  heat  sink-type  wall  liner  with  1 
inch  wall  thickness ,  maximum  wall  temperature  differentials 
are  of  the  order  of  1400°F.  This  results  in  local  stress 
amplitudes  of  240,000  psi.  Since  actual  compressive  stresses 
will  not  exceed  the  yield  point  because  of  thermal  strain, 
cyclic-strain  analysis  will  indicate  the  approximate  number 
of  test  cycles  to  failure.  Cycles  to  failure  are  expressed 
as 


where 


a  = 
a 

a  = 
e 

C  = 
Ra  = 


the  stress  amplitude 
the  endi 
1/2  log 


the  endurance  limit 
100 


100— Rei 

the  area  reduction,  percent 
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APPENDIX  V 

THERMOCOUPLE-PROBE  CORRECTION 


An  analysis  is  made  of  the  radiation  and  conduction  errors 
associated  with  the  platinum  thermocouple  probes  used 
in  the  exit-plane  instrumentation  rake.  Assuming  the  thermo¬ 
couple  bead  to  be  a  hemisphere,  a  summation  of  heat  fluxes 
can  be  written  as 


conv 


^rad 


^oond 


=  0 


(92) 


where 

q  =  convective  heating  flux 
^conv 

q  ,  =  radiation  heat  loss 
rad 

q  ,  =  conductive  heat  loss 
cond 

For  a  hemisphere  with  fluid  properties  taken  at  burner 
design  conditions,  the  Nusselt  number  is  NNu=100.  Taking 
the  surface  area  of  the  bead  as  A=17.5  *  10“6ft2,  the 
convective  flux  can  be  written  as 


q  =hA[T  -T  ]  =  0.193  *  17.5  *  10“6[T.-T  ]  (93) 
^conv  4  w  L  4  w 


where 

T4  =  gas  temperature 
Tw  =  wall  temperature 

The  radiation  heat  flux  can  be  described  by 


q  ,  =  OAfT 
q-ad  w 


where  the  back  radiation  is  ignored  and 


(94) 


<7  =  Stephan-Boltzmann  constant 
(  =  emissivity  =  0.018 

Finally,  the  conduction  is  estimated  by  assuming  the  lead 
wire  to  be  infinitely  long  and  computing  this  heating 
rate  as  a  function  of  time.  From  Reference  11, the  flux 
at  the  interface  between  lead  wire  and  bead  is 
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(95) 


Vtcq] 

qcond  -  JltT?' 

where 

=  wire  temperature  at  infinity  room  temp 
a"  =  thermal  diffusivity  =  k/p 
T  =  time,  sec 

A  =  cross-sectional  area  of  wire 

Table  IX  lists  the  pertinent  characteristics  of  the 
thermocouple  lead  wire. 


TABLE  IX. 

THERMOCOUPLE  LEAD  WIRE 

CHARACTERISTICS 

Tv’-ie 

k=(BTU 

'sec,°R 

-  )  a" 

,ft  (i't2/sec) 

Dia  of 
Wire 
(in. ) 

qcond,BTU/sec) 

Plat. 

10% 

Rhod. 

0.0048 

0.0155 

0.014 

1.99  •  10”7(T  -540) 
w 

(96) 

Plat. 

0.0114 

0.0311 

0.014 

2.98  •  10_7(T  -540) 
w 

(97) 

Solving  the  above  equations  at  time  equal  to  2  seconds 
gives  the  gas  temperature  as  a  function  of  the  junction 
(wall)  temperature  and  is  shown  in  Figure  52. 
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appendix  VI 

TOTAL  TEMPERATURE  DISTRIBUTION  DATA 


This  section  presents  the  experimental  data  associated 
with  the  total  temperature  distributions.  Four  sets  of 
data  are  shown,  corresponding  to  Figures  34  through  37. 

The  value  Tp  represents  the  measured  temperature  and  T„ 
is  the  value  corrected  with  Figure  52. 


TABLE  X. 

DATA  LOG:  ALTITUDE,  100%  POWER 

CONDITION 

Run 

P3 

(psia) 

T3 

(  °R) 

M 

3 

f/a 

Taf 
(  °R> 

t 

tri 

fRl 

T 

R2 

T 

R2 

l 

T 

R3 

T  , 

R3 

e 

(deg  ) 

30 

73 

1075 

0.6 

0.023 

2500 

2300 

2600 

2275 

2590 

2185 

2500 

30 

63 

70 

965 

0.42 

0.022 

2240 

1810 

2000 

1870 

2120 

1810 

2000 

45 

58 

72 

1005 

0.57 

0.022 

2500 

1730 

1900 

1840 

2100 

1765 

2000 

75 

50 

65 

1010 

0.57 

0.022 

2400 

940 

940 

940 

90 

68 

77 

990 

0.45 

0.021 

2300 

1375 

1285 

1285 

225 
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TABLE  XI.  DATA  LOG i  ALTITUDE,  60%  POWER  CONDITION 


Run 

P3 

(psia) 

T3 

(  °R) 

M 

3 

f/a 

af 

(  °R) 

T 

JU 

T 

Rl 

T  „ 

R2 

T 

R2 

T 

_R3 

T 

R3 

0 

(deg  ) 

40 

50 

950 

0.56 

0.021 

2240 

2065 

2340 

2105 

2380 

1920 

2150 

30 

64 

55 

920 

0.57 

0.012 

1700 

1460 

1620 

1540 

45 

59 

53 

980 

0.57 

0.018 

2100 

1650 

1920 

1700 

2000 

1675 

1950 

75 

49 

65 

1010 

0.57 

0.022 

2400 

940 

940 

940 

90 

69 

47 

940 

0.6 

0.02 

2200 

940 

940 

940 

225 

. . - 

TABLE 

XII  . 

DATA  LOG 

:  SEA 

LEVEL, 

100%  POWER 

CONDITION 

P3 

T, 

M 

f/a 

T  4T 

af 

o 

T  , 

R 1 

T 

A’ 

T  „ 

JK2 

T 

_R3 

0 

Run 

(psia) 

(  PI 

3 

(  R) 

T  , 

Rl 

T 

P  2 

T 

R3 

(deg  ) 

43 

168.5 

1157 

0.64 

0.033 

3000 

2770 

3200 

2860 

3300 

2800 

3250 

30 

60 

160 

1050 

0.60 

0.019 

2200 

1785 

2000 

1815 

2050 

1730 

1950 

45 

57 

165 

1040 

0.59 

0.018 

2160 

94  5 

94  5 

945 

75 

48 

161 

1070 

0.59 

0.018 

2200 

1140 

1  100 

1100 

90 

6  5 

167 

1010 

0.54 

.0185 

2200 

1190 

1190 

1190 

225 
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TABLE 

XIII. 

DATA  LOG 

:  SEA 

LEVEL. 

60r/  DOWER 

CONDITION 

T. 

A . 

Pun 

P3 

(psia ) 

T3 

(  °F> 

M 

3 

f/a 

Taf 
(  °R> 

'  tri 
‘ri 

T 

P.2 

T 

R2 

- 1 

T 

R3 

f  , 

R3 

9 

(deg  ) 

44 

140 

1100 

0.66 

0.013 

1900 

1800 

2000 

1835 

2050 

1760 

1950 

30 

61 

140 

1080 

0.60 

0.015 

2000 

174  5 
1950 

1810 

2050 

1745 

1950 

45 

56 

<6 

1085 

0.59 

0.016 

2100 

1285 

1380 

1450 

75 

45 

132 

1000 

0.57 

0.016 

2020 

94  5 

94  5 

945 

90 

67 

133 

1020 

0.67 

0.015 

1950 

940 

940 

940 

225 
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conditions,  i.e.,  pressures  from  1  to  11  atmospheres  and  inlet  air 
temperatures  from  ambient  to  1200°R.  ^  i _ 
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